Draft version February 2, 2008 

Preprint typeset using LAT^^X style cmulateapj v. 21/08/00 



O 
O 

(N 



> 
cn 

m 
o 

o 
6 



THE SPECTRA AND VARIABILITY OF X-RAY SOURCES IN A DEEP Chandra OBSERVATION 

OF THE GALACTIC CENTER 

M. P. MuNo/'^ J. S. Arabadjis/' F. K. Baganoff,'' M. W. Bautz;' W. N. Brandt,^ P. S. 
Broos,^ E. D. Feigelson,'' G. P. Garmire/ M. R. Morris/ and G. R. Kicker^ 

Draft version February 2, 2008 



ABSTRACT 

We examine the X-ray spectra and variability of the sample of X-ray sources with Lx ~ 10'^^ 



10^^ erg 

s~^ identified within the inner 9' of the Galaxy by Muno et al. (2003a). Very few of the sources exhibit 
intra-day or inter- month variations. We find that the spectra of the point sources near the Galactic center 
are very hard between 2-8 keV, even after accounting for absorption. When modeled as power laws the 
median photon index is F = 0.7, while when modeled as thermal plasma we can only obtain lower limits 
to the temperature of kT > 8 keV. The combined spectra of the point sources is similarly hard, with a 
photon index of F = 0.8. Strong line emission is observed from low-ionization. He-like, and H-like Fe, 
both in the average spectra and in the brightest individual sources. The line ratios of the highly-ionized 
Fe in the average spectra are consistent with emission from a plasma in thermal equilibrium. This line 
emission is observed whether average spectra are examined as a function of the count rate from the 
source, or as a function of the hardness ratios of individual sources. This suggests that the hardness of 
the spectra may in fact be to due local absorption that partially-covers the X-ray emitting regions in 
the Galactic center systems. We suggest that most of these sources are intermediate polars, which (1) 
often exhibit hard spectra with prominent Fe lines, (2) rarely exhibit either flares on short time scales or 
changes in their mean X-ray flux on long time scales, and (3) are the most numerous hard X-ray sources 
with comparable luminosities in the Galaxy. 

Subject headings: stars: X-rays — binaries: general — Galaxy: nucleus — novae: cataclysmic variables 
— X-rays: binaries — X-rays: stars 



1. INTRODUCTION 

Recent deep Chandra observations of the inner 9' around 
the super-massive black hole at the Galactic center have re- 
vealed over 2000 individual point-like X-ray sources (Muno 
et al. 2003a). The sources have luminosities between 10"^^ 
and 10'^'^ erg s~^(for a distance of 8 kpc to the Galactic 
center; see McNamara et al. 2000), and thus they proba- 
bly represent some combination of young stellar objects, 
Wolf-Rayet and early O stars, interacting binaries (RS 
CVns), cataclysmic variables (CVs), young pulsars, and 
black holes and neutron stars accreting from binary com- 
panions (low- and high-mass X-ray binaries; LMXBs and 
HMXBs). However, the spectra of the Galactic center 
sources are very hard in the energy range of 2-8 keV. 
Spectra that are similarly hard have only been observed 
previously from magnetically accreting CVs (mCVs) and 
HMXB pulsars. Moreover, seven of the hard sources ex- 
hibit X-ray modulations with periods between 300 s and 
4.5 h, which also suggests that they are magnetized white 
dwarfs or neutron stars (Muno et al. 2003c). These basic 
observations are a good first step toward determining the 
natures of the point sources. However, if their natures can 
be determined conclusively, the large number of sources in 
the field would make it possible to study two important 
pieces of astrophysics: (1) the history of star formation at 
the Galactic center, and (2) the physics of X-ray produc- 
tion in accreting stellar remnants. 



How stars form at the Galactic center is still a mystery, 
because the strong tidal forces and milliCauss magnetic 
fields there should prevent all but the most massive molec- 
ular clouds from collapsing. Nonetheless, it appears that 
star formation has occurred recently, because three mas- 
sive stellar clusters younger than 10^ years old lie within 
~ 30 pc of the Galactic center: IRS 16, the Arches, and 
the Quintuplet (Krabbe et al. 1995; Paumard et al. 2001; 
Figer et al. 1999). However, it is still a matter of debate 
as to whether the star formation is continuous or episodic, 
and whether it occurs only in localized regions or is rela- 
tively uniform throughout the Galactic center. Figer et al. 
(2004) addressed this question by modeling the evolution 
of the population of luminous infrared stars, and concluded 
that the star formation is probably continuous. The X- 
ray sources at the Galactic center could provide an addi- 
tional, independent constraint on the star formation his- 
tory there, because they should be dominated by accreting 
stellar remnants. 

The size of the sample of X-ray sources — an order 
of magnitude larger than the numbers of known LMXBs, 
HMXBs, and magnetic CVs — also makes it a valuable 
database for studying the physics of X-ray production. 
Several outstanding questions could be addressed with the 
current data. If the sample contains large numbers of mag- 
netic CVs, it could be used to determine the duty cycle 
of bright accretion states (e.g., Garnavich & Szkody 
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and the fraction of such systems that exhibit hard spectral 
components (e.g. Ramsay ct al. 2004a). If there is a signif- 
icant number of neutron star HMXBs, it may be possible 
to determine whether material accreted at rates far below 
the Eddington limit can penetrate the neutron star's mag- 
netosphere and reach its surface (Negueruela et al. 2000; 
Campana et al. 2002a; Orlandini et al. 2003). Finally, 
the large sample of Galactic center X-ray sources would 
be useful for identifying systems with unusual properties. 
Previous hard X-ray surveys of the Galactic plane have 
identified several slowly-rotating accreting neutron stars 
and white dwarfs (Kinugasa et al. 1998; Torii et al. 1999; 
Oosterbroek et al. 1999; Sakano et al. 2000; Sugizaki et al. 
2000) , magnetic CVs with extremely strong emission lines 
from He-Ukc Fe (Misaki et al. 1996; Ishida et al. 1998; Tcr- 
ada et al. 1999), and accreting stellar remnants with high 
intrinsic absorption (Patel et al. 2004; Matt & Guainazzi 
2003; Walter et al. 2003). These systems could represent 
resting points for stellar remnants that have not been ob- 
served previously, and are therefore important for calcu- 
lating the formation rate of such remnants in the Galaxy. 

In this paper, we take a further step toward the above 
goals by using the properties of the X-ray emission from 
the point sources near the Galactic center (Muno et al. 
2003a) to constrain better their natures. In Sections 2.1- 
2.3, we examine the spectra of the point sources both indi- 
vidually and averaged together, in order to determine the 
temperatures of the emitting regions. In Section 2.4, we 
search for short-term variability, which is often seen from 
coronal X-ray sources, and long-term variability, which is 
common in some accreting X-ray sources. In Section 3, we 
compare the properties of the observed sources with those 
of known classes of X-ray source. Finally, in Section 4, we 
briefly explore the future prospects for definitively identi- 
fying the natures of these sources. 

2. OBSERVATIONS AND DATA ANALYSIS 

Twelve separate pointings toward the Galactic center 
have been carried out using the Advanced CCD Imaging 
Spectrometer imaging array (ACIS-I) aboard the Chan- 
dra X-ray Observatory (Weisskopf et al. 2002) in order to 
monitor Sgr A* (Table 1). The ACIS-I is a set of four 
1024-by-1024 pixel CCDs, covering a field of view of 17' 
by 17'. When placed on-axis at the focal plane of the 
grazing-incidence X-ray mirrors, the imaging resolution is 
determined primarily by the pixel size of the CCDs, 0'.'492. 
The CCDs also measure the energies of incident photons 
within a calibrated energy band of 0.5 8 keV, with a res- 
olution of 50-300 eV (depending on photon energy and 
distance from the read-out node). The CCD frames are 
read out every 3.2 s, which provides the nominal time res- 
olution of the data. 

The methods we used to create a combined image of 
the field, to identify point sources, and to compute the 
photometry for each source are described in Muno et al. 
(2003a) and Townsley et al. (2003). In brief, for each ob- 
servation we corrected the pulse heights of the events for 
position-dependent charge-transfer inefficiency (Townsley 
et al. 2002b), excluded events that did not pass the stan- 
dard ASCA grade filters and Chandra X-ray center (CXC) 

^ http: //www. astro .psu. edu/xray/docs/TARA/ 

^ http : //tiwa . astro . psu . edu/ users/ chart as/xcontdir/xcont . html 



good-time filters, and removed intervals during which the 
background rate fiares to > 3a above the mean level. The 
final total live time was 626 ks. In order to produce a 
single composite image, we then applied a correction to 
the absolute astrometry of each pointing using three Ty- 
cho sources detected strongly in each Chandra observation 
(compare Baganoff et al. 2003), and re-projected the sky 
coordinates of each event to the tangent plane at the ra- 
dio position of Sgr A*. The image (excluding the first 
half of ObsID 1561, during which the 10^^" erg cm^^ s""'^ 
transient GRS 1741.9—2853 was observed; see Muno et al. 
2003c) was searched for point sources using wavdetect 
(Freeman et al. 2002) in three energy bands: 0.5-8 keV, 
0.5-1.5 keV, and 4-8 keV. We used a significance thresh- 
old of 10~^, which corresponds to the chance probability 
of detecting a spurious source within a beam defined by 
the point spread function (PSF). We detected a total of 
2357 X-ray point sources. Of these, 281 were detected in 
the soft band (124 exclusively in the soft band), and so 
are located in the foreground of the Galactic center. The 
remaining sources, of which 1792 were detected in the full 
band, and 1832 in the hard band (441 exclusively in the 
hard band) are most likely located near or beyond the 
Galactic center. 

We computed photometry for each source in the 0.5- 
8.0 keV band using the acis_extract routine from the 
Tools for X-ray Analysis (TARA).^ We extracted event 
lists for each source for each observation, using a polygo- 
nal region generally chosen to match the contour of 90% 
encircled energy from the PSF, although smaller regions 
were used if two sources were nearby in the field. We used 
a region defined by the PSF for 1.5 keV photons for fore- 
ground sources, and a larger extraction area corresponding 
to the PSF for 4.5 keV photons for Galactic center sources. 
A background event list was extracted for each source from 
a circular region centered on the point source, excluding 
from the event list (i) counts in circles circumscribing the 
95% contour of the PSF around any point sources and 
(ii) the bright, filamentary structures noted by Park et al. 
(2004). The background region was unique for each obser- 
vation. It was chosen to include a fraction of « 1200 total 
counts, where the number of counts from each observation 
was scaled to the fraction of the total exposure time. The 
photometry for the complete sample of sources is listed in 
the electronic version of Table 3 from Muno et al. (2003a). 

We then extracted spectra and background estimates for 
each of the sources from the same regions from which we 
computed the photometry. We summed the source and 
background spectra from all 12 observations. Then, we 
grouped the source spectra between 0.5 8.0 keV so that 
each spectral bin contained at least 20 total counts. Next, 
we computed the effective area function at the position 
of each source for each observation. This was corrected to 
account for the fraction of the PSF enclosed by the extrac- 
tion region and for the time- varying hydrocarbon build-up 
on the detectors.^ We estimated the detector response for 
each source in each observation using position-dependent 
response files that accounted for the corrections we made 
to undo partially the charge-transfer inefficiency (Towns- 
ley et al. 2002a). Finally, to create composite functions 
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Fig. 1. — Example spectra of four bright sources, illustrating the signal-to-noise available from the Chandra observations. The top panels 
have the spectrum in units of detector counts keV"'^ as a function of energy in keV, so that the varying effective area of the detector is 
convolved with the spectrum. The solid histograms represent the best-fit spectra. The bottom panels show the residuals to the fit, in units 
of the difference between the model and the data divided by the uncertainties on the data. Panel a displays a Galactic center source with a 

r = 0.2 power-law spectrum. Panel b illustrates a Galactic center source with a spectrum consistent with a kT = 2 keV plasma. We note 
that this source is far off axis, so its PSF is 50 times larger than for the other sources in the figure. As a result, the relative fraction of 
background counts is larger, and the uncertainty on the net counts also larger. Panel c exhibits a foreground source that can be modeled 
with a two-temperature thermal plasma. Panel d displays a Galactic center source with an Fe line at 6.7 keV with a 2 keV equivalent width, 
as well as lines from ionized Si, S, and Ar. 



of the PSF at 4.5 keV. The spectra were indistinguishable 
for the larger and smaller extraction regions, which con- 
firms that we have successfully removed the background 
emission from the point-source spectra. 

2.1. Spectra of Individual Sources 

We modeled the X-ray spectra of those sources with at 
least 80 net counts, which provided four or more indepen- 

The spectra, response functions, effective area functions, background estimates, and event lists for each source are available from 
http : / /www . astro . psu . edu/users/niel/galcen-xray-data/galcen-xray-data . html. 



for the full data set, we averaged both the response and 
effective area functions, weighted by the number of counts 
detected from each source in each observation. Four ex- 
ample spectra are displayed in Figure l7 

Wc have confirmed that the spectra of the point sources 
were not contaminated by the diffuse X-ray emission in 
the field by repeating the analysis above for a subset of 
sources using an extraction region that enclosed only 50% 
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Fig. 2. — The distribution of the photon indices F, plasma temperatures kT, and absorption columns Nn derived from spectral fits to 
sources with at least 80 net counts. We have only included those sources that are consistent with the spectral models at the 95% confidence 
level. The filled histogram represents the best-fit values. The upper limits from the 90% uncertainty (single parameter of interest) on F and 
ATh are indicated with the dotted histograms, and the lower limit on kT with the dashed histogram. The arrows indicate the numbers of 
sources with best-fit values outside the ranges of the plots. Most of the sources have hard spectra; the median photon index is F = 0.7, while 
we were only able to obtain lower limits on kT from the plasma modes. The median absorption column is 6 x 10^^ cm~^. 



dent spectral bins. To provide a rough characterization 
of the spectrum, we used either a power-law or thermal 
plasma continuum absorbed at low energies by gas and 
dust. To model the thermal plasma, we used mekal in 
XSPEC (Mcwc, Lcmcn, & van den Oord 1986). Wc assumed 
that the elemental abundances were 0.5 solar, which is con- 
sistent with the values derived from the average spectra of 
the point sources (see Section 2.3.2), as well as the Fe 
abundances often observed from CVs (e.g., Done & Os- 
borne 1997; Fujimoto & Ishida 1997; Ishida ct al. 1997).^ 
We accounted for gas absorption using the model phabs, 
and the dust scattering using a modified version of the 
model dust in which we removed the assumption that the 
dust was optically thin. The column depth of dust was set 
to T 0.485 • A^H/(102^cm-2)^ and the halo size to 100 
times the PSF size (Baganoff et al. 2003). In Table 2 we 
list the parameters of the best-fit spectral models: the col- 
umn densities A^h. cither the power-law slope F or the tem- 
perature fcT, the observed and de-absorbed 2-8 keV fluxes, 
and the reduced ■ The uncertainties are 90% confidence 
intervals (Ax^ = 2.71). We also indicate sources from 
which the spectra should be viewed with caution, of which 
there are three categories: confused sources for which the 
radii of their PSF overlap those of a nearby source by more 
than 25%, sources that were near chip edges, and sources 
with variability (Section 2.4). About 25% of the sources 



are flagged in this manner. 

We consider a spcictrum to be adequately reproduced 
by a model if the chance probability of obtaining the de- 
rived value of is greater than 5%. Of the 566 sources 
that we modeled in this manner, the spectra of 470 could 
be modeled with an absorbed power law, and 469 could 
be modeled with an absorbed, coUisionally-ionized plasma. 
Both spectral models were consistent with the data for 440 
sources, because of limited statistics and the small band- 
pass over which photons arc detected (for most sources. 
Galactic absorption prevents photons < 2 keV from reach- 
ing the detector, while the effective area of the ACIS-I is 
small above 8 keV). Only 30 sources could be modeled with 
a power law but not a thermal plasma, while 29 sources 
could be modeled by a thermal plasma but not a power 
law. Finally, 67 spectra deviated from both continuum 
models at the 95% level. About 25 of these resemble sta- 
tistical fluctuations. 30 are bright sources that appear to 
require multiple continuum components, and 12 exhibited 
strong iron emission between 6 and 7 keV. 

For individual sources, the uncertainties on the spec- 
tral parameters are rather large. However, it is possible 
to draw some general conclusions about the population 
of X-ray sources from Table 2. In Figure 2, we plot the 
distributions of the best-fit absorption columns and pho- 
ton indices or temperatures for all of those sources that 



® We note that the abundance parameter merely measures the relative strengths of the lines and the continuum. If the continuum is non-thermal 
or the lines are produced by photo-ionization, the abundance parameter will not measure the physical abundances of metals in the plasma. 
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Fig. 3. — Plot of the best-fit power-law spectral index (F) as a 
function of the observed 2-8 keV flux. Representative uncertainties 
on F arc indicated for a few points. There is no correlation between 
F and the flux. 
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Fig. 4. Comparisons of the chance probability that an iron line 
would be detected from sources, derived using Markov-Chain Monte 
Carlo simulations and using the F-distribution. The F-distribution 
generally over-estimates the chance probability of the line, which 
would cause the significance of a line to be under-estimated. 




were adequately fit with the simple absorbed continuum 
model. As was mentioned in Muno et al. (2003a), 265 out 
of 470 of the point sources have power-law spectra with 
r < 1, even after accounting for the absorption column. 
However, only 77 sources have 90% confidence limits for 
which r < 1, and some of the apparent hardness of their 
spectra is due to line emission from Fe between 6 7 kcV 
(see Section 2.2). Not surprisingly, the hard sources also 
have high best-fit temperatures under the thermal plasma 
model. However, because of the poor statistics and small 
bandpass, the temperatures are unconstrained in one-third 
of the sources, and there are only 21 sources with 90% 
lower limits on the temperatures that are > 10 keV. The 
slopes of the spectra are not correlated with the intensities 
of the sources (Figure 3). 

The median absorption column toward the the sources 
is 6 X 10^^ cm~^ for a power-law continuum, and 11 x 
10^^ cm~^ for a thermal plasma continuum. The median 
value for the power-law continuum is identical to the col- 
umn that is inferred from the K-hand extinction toward 
Sgr A* (see Tan & Draine 2003, for a summary), while that 
for a thermal plasma continuum is twice as high. The dif- 
ference between the median values results from the facts 
that (1) for a faint, hard source the values of Ah and the 
spectral slope cannot be determined independently, and 
(2) the thermal plasma continuum, can appear no harder 
than a r « 1.5 power law. Therefore, the thermal plasma 
models produce a higher median value of A^h- Under ei- 
ther model, 30% of the sources are inferred to have a ab- 
sorption column of > 10^"^ cm~^. The one source that 
appears to be absorbed by more than 10^'* cm~^, CXOGC 
J174539. 3-290027, is relatively faint, and the spectrum is 
poorly constrained. 

Finally, we note that these spectra have been modeled 
assuming that all of the X-rays from a given source are 
absorbed by a single column of material. If we assume 
instead that a fraction of the X-ray emitting region is ab- 
sorbed by a higher column (so-called "partial covering" 
models) an acceptable fit can be obtained for an arbi- 
trary range of continuum shapes, because the bandpass 
over which we measure the spectrum is limited. 



2.2. Iron Emission 

Visual inspection of the 67 sources for which the simple 
continuum failed to reproduce their spectra indicates that 
« 20% exhibit residuals between 6-7 keV that may repre- 
sent line emission from Fe. Therefore, we have performed 
a uniform search for Fe line emission from the brightest 
sources with hard spectra. We selected only those sources 
with more than 160 net counts, because we found that 
fainter sources were unlikely to provide more than one 
spectral bin between 6-7 keV. Likewise, we selected only 
those sources that were best-fit by absorbed power laws 
with r < 5, because sources with steeper spectra were 
typically background-dominated above 6 keV. There were 
183 sources that met both of these criteria. We modeled 
each source with an absorbed power law plus a Gaussian 
line that was fixed at either 6.4 keV to search for low- 
ionization ("neutral") Fe emission, or at 6.7 kcV to search 
for He-like Fe. The widths of the lines were fixed at 100 
eV, to account for the fact that both the low-ionization 
and He-like lines are actually a blend of multiple transi- 
tions (e.g., Nagasc ct al. 1994). 

To evaluate the significance of the added line, we com- 
puted a statistic / from the reduction in provided by 
the more complex model: 

(1) 

where x? and Xs are the values with and without the 
line, and Vc and Vs are the numbers of degrees of freedom 
for the fits {vs — — 1, here). Unfortunately, because 
the null result of our more complex model, a line of zero 
flux, lies on a boundary of the parameter space we are 
considering (i.e., an emission line with negative flux is 
not physical, so the hard lower limit to the line flux is 
0), / is not distributed according to the distribution 
(Protassov et al. 2002). Therefore, we have simulated 
the expected distribution of / for each of our sources us- 
ing the Markov-Chain Monte Carlo technique described 
in Arabadjis et al. (2003). This technique allows us to 
evaluate the chance probability of observing a value of 
/ given (1) the statistical distributions of A^h, F, and 
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FlG. 5. — The strength of line emission at left panels: 6.4 keV from low-ionization Fe, and right panels: 6.7 keV from He-like Fe. The top 
panels display histograms of the equivalent widths of the iron lines. The rightward-pointing arrow denotes the the four sources that may have 
iron ciuissioii with equivalent width > 2500 keV. Most of these sources have very steep (F > 4.5) spectra, and so the excess emission between 
6—7 keV could represent either a hard continuum component, or poor background subtraction. However, the line from one of the sources with 
6.4 keV iron emission with equivalent width > 2500 cV appears to be real (CXOGC J174617. 2— 285449). The bottom panels illustrate how 
the measured equivalent widths are correlated with the photon flux from a source. Iron emission can only be detected in the fainter sources 
if it has a high equivalent width. 
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Fig. 6.- - Combined spectrum of all Galactic center point sources 
with fewer than 500 net counts, before background subtraction. The 
background spectrum is taJsen to be that of the diffuse emission, and 
is indicated with the dashed curve. 



the normalization from the fit, (2) the method used to 
group the spectral bins, and (3) the spectrum assumed for 
the background subtraction. Between 100 1000 simula- 
tions^ were computed for each source. We present in Fig- 
ure 4 a comparison of the chance probabilities derived from 
the theoretical i^-distribution and from the Monte Carlo 
simulations. In general, the theoretical i^-distribution sig- 
nificantly under-estimates the significance of a line (the 
chance probability is deemed to be too large), although in 
10% of the cases the theoretical distribution over-estimates 
the significance. This demonstrates the necessity of per- 
forming simulations in order to estimate the significance 
of an added line feature. 

We have listed those sources with significant line emis- 
sion in Table 3. We consider a line to be detected if it 
has less than a 1% chance probability of arising from ran- 
dom variations in the continuum, and if the fit produced 
xl/i^c < 2. In total, 35% (64 out of 181) of the sources that 
we examined have significant Fe line emission. We find 
that 25% of sources exhibit a significant line near 6.7 keV, 
while 16% of sources exhibit emission at 6.4 keV. There are 
15 sources with significant lines at both 6.4 and 6.7 keV, 
but these tend to be faint sources in which it is not possi- 
ble to constrain the line energy. We plot histograms of the 
equivalent widths and the equivalent widths as a function 
of the flux in Figure 5. For both species, the equivalent 
widths range from 200 eV to 5 keV. The median equivalent 
widths were 370 eV for the 6.4 keV Fe line, and 530 eV 
for the 6.7 keV Fe line. Six sources have lines with equiv- 
alent widths greater than 1 keV that appear upon visual 
inspection to be real. There are two systems with ap- 
parent equivalent widths > 10 keV, but these have steep 
continuum spectra with F > 4.5, and the excess emis- 
sion between 6-7 keV represents either a hard continuum 
component or poor background subtraction. When line 
emission is not detected, the median 1-a upper limits to 
the line equivalent widths are 220 eV for He-like Fe, and 
240 cV for low-ionization Fe. 




Fig. 7.- Combined, background-subtracted spectrum of the 
Galactic center point sources with fewer than 500 net counts. The 
best-fit model of a power-law plus 10 Gaussian lines is over-plotted. 
The expected centroid energies of lines from Si, S, Ar, Ca, and Fe 
are indicated. 



2.3. Combined Spectra of Point Sources 

In order to understand the average spectra of the Calac- 
tic center point sources, wc summed the spectra of sub- 
groups of the individual point sources. We selected only 
those sources that were not detected below 1.5 keV with 
wavdetect, as these are most likely to lie near the Galac- 
tic center. Wc also excluded sources brighter than 500 net 
counts, because these sources provided individual spectra 
of good quality. We computed average effective area and 
response functions by averaging those from the individ- 
ual sources, weighted by the number of counts from each 
source. We estimated the average background by extract- 
ing a spectrum from the rectangular region that traced the 
orientation of the ACIS-I detector during the 500 ks series 
of observations from 2002 May to June."'^'^ We excluded 
from the background spectrum events that fell within cir- 
cles circumscribing the 95% contour of the PSF around 
any point sources. 

In Figure 6, we display the summed spectra of the Galac- 
tic center point sources with fewer than 500 net counts. 
The spectrum below 2 keV is dominated by the diffuse 
emission from the Galactic center. In Figure 7 we display 
the background-subtracted spectrum. The instrumental 
Ni line at 7.5 keV is absent from the spectrum, which in- 
dicates that the background subtraction was successful. 
Lines of Si, S, and Ar are also weak or absent in the spec- 
trum of the point sources, in contrast to that of the diffuse 
emission (dashed line in Figure 6). Prominent lines from 
He- like and H-line Fe are evident at 6.7 and 6.9 keV, while 
weak, fluorescent K-a emission from low-ionization Fe is 
evident at 6.4 keV. 

We modeled this spectrum using two approaches. First, 
we modeled the emission phenomenologically, using a 
power-law continuum and Gaussian line emission from Si, 



^ We first computed 100 simulations. For many of the sources, the / values from the simulations exceeded the observed / value more than 
once, indicating the significance of the line was 99%. For the rest, we ran 1000 simulations, to establish whether the line was at least 99% 
significant (i.e., < 10 / values exceeding the observed one). 

We did not simply sum the background spectra obtained for individual sources, because doing so would have double-counted events from 
background regions that overlapped. 
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Fig. 8. — Ratio of the spectra of point sources with between 80 
and 500 net counts to those with < 80 net counts. The ratio is fiat 
between 2.5—6 IceV, which indicates that the continuum shapes are 
nearly identical. The deviations at 2.4 keV are due to stronger S 
emission in the brighter sources. The deviations between 6-7 keV 
are due to the range of strengths in the neutral, He-like, and H-like 
lines of Fe. 

S, Ar, Ca, and Fc. The lines we included were chosen 
by examining whether lines from the strongest transitions 
from Table 1 in Mewe, Gronenschild, & van den Oord 
(1985) significantly improved the residuals when compar- 
ing the model to the data. For the final model, we placed 
lines at the energies expected for the He-like n = 2 — 1 
transitions of Si, S, Ar, Ca, and Fe; the He-like n = 3 — 1 
transitions of Si and S; the H-like n = 2—1 transitions of Si, 
S, Ar, and Fe; and low-ionization Fe K-a at 6.4 keV. This 
model allowed us to measure and compare the eqtiivalent 
widths of the lines. We also used a model consisting of two 
thermal plasma components, each of which was absorbed 
by a separate column of gas. This model is identical to the 
model used for the diffuse emission by Muno et al. (2004). 
We also note that this model is qualitatively similar to the 
multi-temperature, multi-absorber models typically used 
to model the accretion shocks in magnetized CVs (e.g. 
Ramsay et al. 2004a). 

Several assumptions were required for the models to re- 
produce the data. First, we only applied the model be- 
tween 1.0 — 8.0 keV. Below this energy range the pho- 
ton counts are dominated by foreground diffuse emission, 
while above this range the ACIS-I has a small effective 
area. Second, when modeling individual lines with Gaus- 
sians, the widths of the lines from He- like Si, S, and Fe and 
low-ionization Fe were allowed to be as large as ~ 70 eV 
to account for the fact that the lines are blends of sev- 
eral transitions that cannot be resolved with ACIS. Third, 
we allowed for a < 1% shift in the energy scale in each 
spectrum because of uncertainties in gain calibration of 
our CTI-corrected data. When fitting Gaussians to the 
He-like transitions and the 6.4 keV line of Fe, the line cen- 
troids were varied one-by-one until they achieved best-fit 
values, and then frozen. When using plasma models, the 
red-shift parameter was used to change the energy scale 
in a similar manner. Finally, a 3% systematic uncertainty 
was added in quadrature to the statistical uncertainty in 
order to account for uncertainties in the ACIS effective 
area. 

To account for absorption in each model we assumed (1) 




3 4 5 6 
Energy (keV) 



Fig. 9. — Combined spectrum of Galactic center point sources 
with fewer than 500 net counts and sorted according to their hard 
color. Although the observed shape of the continuum varies strongly 
with hardness ratio, lines from He-like and H-like Fe are present in 
all three groups of sources. The ratios of the Fe lines axe consistent 
with a kT » 8 keV plasma in all cases. This suggests that varia- 
tions in the local absorption column toward the sources cause the 
variations in the apparent hardness of the Galactic center sources. 



that the entire region was affected by one column of mate- 
rial that represents the average Galactic absorption (mod- 
eled with phabs in XSPEC) and (2) that a fraction of each 
region was affected by a second column that represents ab- 
sorbing material that only partially-covers the X-ray emit- 
ting region (modeled with pcf abs). This partial-covering 
absorption model produces a low-energy cut-off that is 
less steep than that which would be produced by a sin- 
gle absorber. The model can roughly account for the fact 
that both the point sources and absorbing material are dis- 
tributed along the line of sight. The mathematical form 
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of the model was 

where a{E) is the energy-dependent absorption cross- 
section, A^H is the absorption column, -/Vpc,H is the partial- 
covering column, and /pc is the partial-covering fraction. 
Dust scattering was not included, because when modeling 
the spectrum its optical depth was degenerate with the 
partial-covering fraction /pc- 

2.3.1. Phenomenological Model 

Since the natures of these point sources are uncertain, 
the most straightforward way of modeling their spectrum 
is with an absorbed power-law continuum and Gaussian 
line emission. The average spectrum of soiirccs with fewer 
than 500 net counts is displayed along with the model spec- 
trum in Figure 7. The model parameters are listed in the 
first column of Table 4. After some initial tests, we found 
that /pc was poorly constrained, but the best-fit values 
were near 0.95. Wc therefore fixed /pc to this value. The 
remaining parameters were allowed to vary. The total ab- 
sorption column. Nil + -^pc.H » 9 x 10^^ cm~^, is slightly 
higher than the expected Galactic value. Simulations in 
XSPEC indicate that this is probably because we did not 
include dust scattering, which produces about 30% of the 
total absorption. The inferred continuum is flat, with pho- 
ton index T = 0.8, which is similar to the median value 
from the individual sources, T = 0.7. Finally, the equiva- 
lent width of the He-like Fe line is 400 eV, which is sim- 
ilar to that observed from individual bright sources. The 
strength of the neutral Fe line is considerably lower than 
those detected from the individual sources, although this is 
not surprising since fewer of the individual sources exhibit 
6.4 keV lines than do 6.7 keV lines (see Section 2.2). 

We then examined how the average spectra of the point 
sources varied with flux. We compared the summed spec- 
tra of sources with fewer than 80 net counts to those 
with 80-500 net counts, because these two groups of 
sources produce nearly the same numbers of net counts 
(» 7.5 X 10^). The best-fit parameters of our phenomeno- 
logical model of these spectra are listed in the second two 
columns of Table 4. Both the absorption column and pho- 
ton index were slightly larger in the bright sources. How- 
ever, the continuum shape looks very similar by eye. To 
highlight this, in Figure 8 we plot the ratio of the aver- 
aged spectra of sources with 80 500 net counts to that of 
sources with < 80 net counts. The only differences in the 
continuum spectra are above 7 keV, which may indicate 
that the faint sources produce slightly more high-energy 
flux. Differences in the equivalent widths of the line emis- 
sion are also evident in Figure 8. The equivalent width of 
the Fe XXV He-a emission was 30% higher in the fainter 
sources {3.7a), and the equivalent width of neutral Fe K-a 
was a factor of 2 lower in the faint sources (4.9cr). There is 
also a factor of 3 less S He-a emission in the faint sources 
(4.3(t). We also split the data into smaller flux intervals, 
and found similar results. Thus, the continuum spectra 
of the point sources change very little with intensity, but 
the S and Fe lines range over a factor of > 2 in equivalent 
width. 

We also examined the combined spectra as a func- 
tion of the hard color, which provides a measure of the 
steepness of the spectrum. The hard color is defined as 



HR = {h — s)/{h + s), where h is the number of counts 
in a hard band from 4.7-8.0 keV, and s is the number 
of counts in a softer band from 3.3-4.7 keV. We divided 
the data using three ranges in hard color, guided by the 
expected power-law index from simulations with PIMMS 
(Muno et al. 2003a): HR > 0.2 for sources with F f« 
spectra, —0.1 < HR < 0.2 for sources with F 1.5 spec- 
tra, and HR < —0.1 for sources with F > 2 spectra. The 
spectra are plotted in Figure 9, and the best-fit parame- 
ters for the phenomenological model are listed in the last 
three columns of Table 4. We find that the sources with 
larger HR have continuum spectra that are intrinsically 
fiatter, at least given our model for the interstellar and in- 
trinsic absorption. We also find significant changes in the 
line strengths, although there is no monotonic trend with 
HR. 

2.3.2. Two-kT Plasma Model 

We next modeled the spectra of the point sources as 
originating from two thermal plasmas, each of which was 
absorbed by a separate component as parameterized in 
Equation 2, as well as emission from low-ionization Fe at 
6.4 keV. The parameters of the best fit models to the aver- 
age spectra for the groups of sources used in the previous 
section are all listed in Table 5. The values of xl are gen- 
erally larger under the plasma models than under the phe- 
nomenological models, because the plasma models predict 
too little flux above 7 keV. However, the plasma models 
do provide a good qualitative description of the data. The 
weak lines from Si, S, and Ar require cool plasmas with 
kTi « 0.5 keV; this component appears to be visible only 
in sources with HR < —0.1. The abundances of Si and 
S appear to be significantly below solar values in the soft 
sources, where the lines are most clearly detected. The 
He-like and H-like Fe lines require hotter plasmas with 
kT2 « 8 keV. The abundances of Fe are also generally 
about 50% solar, although they are consistent with the 
solar ones in the faintest and softest sources. 

The temperatures of these plasma components are very 
similar to those inferred from the diffuse emission, because 
emission lines from the same set of ions are present in both 
spectra (Muno et al. 2004). However, in contrast to the 
diffuse emission, the cooler plasma in the point sources 
is heavily absorbed and contributes little to the observed 
flux. 

More importantly, the presence of prominent He-like 
and H-like Fe emission from kT = 7 — 9 keV plasma in 

all groups of sources suggests that the X-ray emission is 
produced by similar physical mechanisms. 

2.4. Search for Variability 

We searched for variability in the 0.5-8.0 keV bandpass 
from the entire sample of X-ray sources from Muno et al. 
(2003a), in order to identify flux changes that occurred 
within single observations, and long-term variations in the 
mean flux between observations. 

2.4.1. Short-term Variability 

To search for short-term variability, we applied a 
Kolmogorov-Smirnov (KS) test to the un-binned arrival 
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Fig. 10. — Example light curves of variable sources in the field. To compute the resulting photon fluxes, the total counts received from each 
source as a function of time have been divided by the mean value of the detector effective area from each observation. 
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Fig. 11. — Summary of the amplitudes and durations of ttic variability in sources detected toward the Galactic center. The top left panel 
displays a histogram of the ratio of the maximum to the minimum flux from each variable source. Sources that were detectable at minimum 
are indicated with solid histogram, while lower limits on the variability amplitude are illustrated with the open histogram. The bottom left 
displays the variability amplitude (or lower limits thereto) as a function of the mean photon flux from the source. Not surprisingly, small- 
amplitude variability can only be detected if a source is bright. The top right panel displays a histogram of the durations of flare-like events. 
These range from < 1 ks to nearly 100 ks. The bottom right panel illustrates that the observed flare duration is not strongly correlated with 
the mean photon flux from a source. 



times of the events during each observation. Before per- 
forming the search, we removed events flagged as poten- 
tial cosmic ray after-glows. We also excluded data received 
near the edges of the detector chips, and data from the first 
part of ObsID 1561 for sources that were within 5.5' of the 
bright transient GRS 1741.9-2853 (Muno et al. 2003b). 
If the cumulative distribution of the arrival times differed 
from a uniform distribution (which would imply a constant 
flux) with greater than 99.9% confidence in any observa- 
tion, we considered the source to vary on short time-scales. 
We find that 18 foreground sources and 21 Galactic cen- 
ter sources are variable on short time scales according to 
the KS test. We list these sources in Table 6. Examples 
of short-term variability are shown in the bottom three 
panels of Figure 10. 

To characterize the duration and amplitude of the vari- 
ability, we have applied the"Bayesian Blocks" algorithm 
of Scargle (1998) (see also Eckart et al. 2004). The al- 
gorithm is based on a parametric maximum-likelihood 
model of a Poisson process that divides the data into se- 
quential segments, each of which has a constant count 
rate. The segments were identified by dividing the events 
into sub-intervals, and computing the odds ratio that 
the count rate has varied. If variability was found, then 
each interval was split further into sub-intervals, in order 
to track the structure of the variability. We found that 
by using an odds ratio corresponding to a 67% chance 



that the variation is real, we could identify changes in 
the flux from all but four of the variable sources. If the 
Bayesian Blocks code identified only two intervals with 
differing count rates, then the variability was classified as 
a "step" function. If it identified more than two inter- 
vals with differing count rates, we defined the variability 
as a flare; a search of the data revealed no instances of 
dips, which are unlikely to be detected in sources with 
low count rates. We computed the background-subtracted 
maximum and minimum count rates for each variable ob- 
servation, and divided them by the mean value of the 
effective area function for that source to convert them 
into photon fluxes. The maximum and minimum fluxes 
are listed in Table 6. We also list the durations of the 
flares in kiloseconds, and either the ratio between the 
maximum and minimum flux, or the lower limit thereto if 
the baseline flux is an upper limit. The photon fluxes in 
the table can be converted to energy fluxes by the factors 
1 ph cm-2 s-i== 3 X 10"^ erg cm'^ s"! (0.5-8.0 keV) 
for foreground sources, and 1 ph cm^^ s^^= 8 x 10~^ 
erg cm~^ s~^ (2-8 keV) for Galactic center sources. The 
peak luminosities of the variable Galactic center sources 
typically range from 6 x 10^^ to 1.7 x 10^^ erg s~^. How- 
ever, one flare from CXOGC J174552. 2- 290744 consists of 
5 photons received in 100 s and has a peak luminosity of 
erg s ^ None of the events in the apparent flare were 
flagged as cosmic rays by the CXC pipeline, and tests of 
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curs with count rate r/j lasting th, then the total number 
of counts in each interval follows the Poisson distribution. 
Therefore, the joint probability that the measured base- 
line count rate Ni / ti is less than the measured flare count 
rate Nh/th is 

P{Nh > Nith/ti) = 
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This probability represents the chance that a flare of am- 
plitude Vh/ri would be detected. 

The median net counts from the sources in the cata- 
log of Muno et al. (2003a) was 49, with a background 
of 52 counts. These values translate to count rates of 
n = 7.8 X 10~^ count s"^ net and 6 = 8.3 x 10~^ count s~^ 
background. If we use n = n -|- 6 in Equation 3, a 36 ks 
flare during the long 150 ks observation coifld be detected 
with an amplitude a factor of « 10. Such a flare could be 
detected from half of the sources in our sample. On the 
other hand, a flare that reaches twice the quiescent flux 
level for 36 ks could only be detected if the quiescent count 
rate was r; = 2 x lO"'' count s"'^. Only 17 sources are this 
bright, so such a small- amplitude flare would generally be 
unobservable. Not surprisingly, all of the short-time scale 
variability with amplitudes < 3 are long duration, step-like 
changes in flux. 

2.4.2. Long-term Variability 

To search for long-term variability, we computed the 

value of for the photon fluxes in each observation un- 
der the assumption that the mean rate was constant. We 
computed the approximate total (source plus background) 
photon flux by dividing the total number of counts de- 
tected by the live time and the mean value of the eff'ec- 
tive area function. To compute a net flux, we then sub- 
tracted a background count rate, which was estimated in 
the same manner as for the spectrum. We considered a 
source as variable if the photon fluxes both before and af- 
ter background subtraction were inconsistent with a con- 
stant mean value with more than 99% confidence. ^-'^ We 
excluded sources with short-term variability when search- 
ing for long-term variability. We also excluded data from 
the first part of ObsID 1561 for sources that were within 
5.5' of the bright transient GRS 1741.9-2853. Long-term 
variability is illustrated in the top three panels of Fig- 
ure 10. 

We find that 20 foreground sources and 77 Galactic cen- 
ter sources vary on long time scales. We list in Table 7 
the minimum and maximtim background-subtracted pho- 
ton fluxes for the variable soiirc:es. We present a histogram 
of the ratio of the maximum to minimum fluxes in the top 
panel of Figure 12. Most of the ratios are upper limits, be- 
cause the sources are not detected at their minimum flux 
levels. The bottom panel illustrates the variability ampli- 
tudes as a function of the mean intensity of each source. 
There is no apparent correlation between the amplitude 
and intensity of the source. 

We can quantify our sensitivity to long-term variations 
in the same manner as for short-term variability, using 

11 The requirement on the total flux was designed to ensure that the variability was not due systematic changes in the background estimate. 
Such changes occurred where there were gradients in the diffuse emission, because the regions in which the background were estimated were 
not identical for each observation (Section 2). 



Fig. 12. — Summary of the amplitudes and durations of the vari- 
ability observed between observations in sources detected toward the 
Galactic center. The top panel displays a histogram of the ratio of 
the maximum to the minimum flux from each source with long-term 
variability. Sources that were detectable at minimum are indicated 
with solid histogram, while lower limits on the variability amplitude 
are illustrated with the open histogram. The bottom displays the 
variability amplitude (or lower limits thereto) as a function of the 
mean photon flux from the source. Once again, small-amplitude 
variability can only be detected if a source is bright. 

events from elsewhere on the detector that were flagged as 
cosmic rays indicate that they deposit their energy on time 
scales of < 20 s, so we consider the flare real. No source 
exhibits flares similar to those seen about once a day from 
Sgr A*, with durations of 1 h and Lx ^ lO''^ erg s~^. 

Figure 11 displays histograms of the variability ampli- 
tude and the durations of the flares. Nearly half of the 
variability has a peak flux over 10 times the quiescent 
level. In the bottom-left panel, we plot the amplitude of 
the variability as a function of the mean flux from each 
source. The amplitudes of the variability are not a strong 
function of the mean flux from the source. However, we 
are unable to detect the faintest sources when they are 
in their low-flux states, so for these sources we only can 
report lower limits to the variability amplitude. The flare 
durations are spread fairly evenly, with a median duration 
of about 20 ks. As can be seen from the bottom- right 
panel of Figure 1 1 , the flare durations show no correlation 
with the mean flux from a source. 

In order to quantify our sensitivity to short-term vari- 
ations, we need to examine the probability that a change 
in count rate could be detected. If we assume a baseline 
count rate ri persists for a time ti, and that a flare oc- 
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Equation 3. The most extreme form of long-term variabil- 
ity is that of a source that is bright for some portion of the 
observations lasting a total time t^, and decreases below 
the background level for the remaining observations lasting 
a time ti. We therefore assume that the total coimts from 
a source is consistent with the background 6 = 8.3 x 10~^ 
count s~^ during the time ti when it is faint. If the source 
is "off" during one of the 12 ks observations, we could 
detect this decrease in count rate with 99% confidence if 
during the remaining 614 ks of observations the soiu c;e was 
brighter than 3 x 10"** net count s~^. Approximately 9% 
of the sources are this bright. We could detect variabil- 
ity from a source that is "off" during all but the 500 ks 
monitoring campaign (2002 May-June) if the count rate 
at maximum was 8 x 10~^ count s~^, which is valid for 
half of the sources searched. 

3. DISCUSSION 

The large number of sources detected toward the Galac- 
tic center is most likely a product of the large density of 
stars there. The 17' by 17' field spans a physical distance 
of 20 pc in projection from Sgr A*, and therefore probes 
the inner regions of the Nuclear Bulge that was studied 
extensively by Launhardt, Zylka, & Mezger (2002). Using 
their models, we estimate that 1.3 x 10^ Mq of stars lie 
within a cylinder of radius 20 pc and depth 440 pc that 
is centered on the Galactic center. Thus, our observation 
encompasses up to 0.1% of the total Galactic stellar mass, 
which is ^ 10"'^"'^ Mq. However, we have also found that 
the surface density of the X-ray sources falls off as 
away from Sgr A* , so it is possible that most of the X-ray 
sources lie in an isothermal sphere of radius 20 pc (Muno 
et al. 2003a). Such a sphere would contain 3 x lO'' Mq of 
stars, or 0.03% of the Galactic stellar mass. For compari- 
son, the shallower survey carried out by Wang et al. (2002) 
covered the entire Nuclear Bulge (albeit with a factor of 
5 less sensitivity), and thus sampled 1% of the mass of 
stars in the Galaxy. The stellar density at the location of 
the X-ray sources is between 240-900 Mqpc"^ (for a 20 
by 440 pc cylinder and a 20 pc sphere, respectively), com- 
pared to 0.1 Mq pc~^ in the local stellar neighborhood 
(Binney & Tremaine 1994, p. 16). We will keep these 
numbers in mind as we consider the likely natures of the 
Galactic center point sources. 

The sample identified as part of the Chandra obser- 
vations of the Sgr A* field is unique, because the long 
exposure time allows us to detect faint sources (Fx = 
(3 - 100) x 10^15 erg cm'^ s'^, 2-8 keV), whereas the 
strong diffuse X-ray emission and the high absorption to- 
ward the Galactic center prevent us from observing X-ray 
sources unless they are prominent in the 4-8 keV band. 
Prior to Chandra, the most sensitive hard X-ray survey 
of the Galactic plane was taken with ASCA. That sur- 
vey identified only 163 sources, with a detection limit of 
« 3 X 10-13 erg cm'^ s"! (2-10 keV; Sugizaki et al. 2001). 
The Galactic center sources are on average much harder 
than those detected in the ASCA survey. The brighter 
ASCA sources had a median photon index of F = 2.5, 

As we discuss in Section 2.3, if there is local X-ray absorption that affects only a fraction of the emitting region, the inferred spectrum 
can seem artificially fiat, and the inferred luminosity would be too low. Therefore, the hardness ratios and photon fluxes could potentially 
be misleading. This is nearly impossible to avoid. Indeed, the spectral models we applied to the individual sources suffer from the same 
shortcoming, because we also assume that the emitting region is absorbed uniformly. Using hardness ratios and photon fluxes is the best 
option, since unlike the parameters of the spectral models, the former can be derived for almost all of the sources in our sample. 



with only 15% of the sources having F < 1, while the 
fainter Chandra sources have a median F = 0.7 (Figure 2). 
The difference in hardness of the two samples is probably 
a selection effect caused by the high absorption and the 
strong diffuse emission toward the Galactic center. The 
X-ray sources that we have studied in this paper proba- 
bly sample only the hardest examples of the population 
identified with ASCA. 

Likewise, Chandra observations of globular clusters have 
identified a couple hundred X-ray sources with Lx = 
1029 _ 1033 erg s-i(Grindlay et al. 2001; Pooley et al. 
2002; Becker et al. 2003; Heinke et al. 2003a,b). These lu- 
minosities overlap those inferred for the sources near the 
Galactic center in our sample. However, only a few of the 
~ 60 globular cluster sources with spectral information are 
best modeled by a F < 1 power-law. Most have steeper 
spectra that are consistent with F > 1.5 power-laws or 
kT w 1 — 20 keV thermal plasmas. 

In addition to the hardness of their spectra, the X-ray 
sources detected toward the Galactic center share several 
other interesting properties. Most notable is line emission 
from low-ionization. He-like, and H-like Fe (Figures 5 and 
7). On average, the low-ionization Fc lines have equiva- 
lent widths of 100-230 eV, while the He-like Fe lines have 
equivalent widths of 350-450 eV. The strengths of these 
lines range over a factor of two when considering sources 
with a range of intensity and spectral hardness (Table 4) . 
However, in all cases the average ratios of the He-like and 
H-like lines are consistent with those expected from a ther- 
mal plasma of fcT « 8 keV (Table 5). The presence of these 
Fe lines in a large fraction of the sources suggests that they 
could be dominated by a single population of sources. 

However, the emission from such a plasma should pro- 
duce a much steeper continuum spectrum, with F « 1.5 
instead of F » 0.7. Unfortunately, it is not possible to 
determine unambiguously the physical process producing 
the X-ray emission from the continuum and iron lines 
alone. For instance, if the X-ray emitting regions are par- 
tially absorbed by material local to the X-ray sources, the 
observed spectra can be much harder than the intrinsic 
ones. Alternatively, the line emission could be produced 
in photo-ionized plasmas, although the large equivalent 
widths of the lines indicates that the continuum emission 
exciting them must not be observed directly (e.g. Mukai 
et al. 2003). In either of these cases, the intrinsic X-ray 
luminosity woidd be significantly higher than is inferred 
from the flux received with Chandra. 

While the average spectral properties provide an 
overview of the characteristics of the X-ray sources near 
the Galactic center, it is still important to examine the 
properties of individual sources to determine how vari- 
ous classes of sources contribute to the population there. 
Therefore, in Figure 13 we display the spectra and inten- 
sities of individual sources by plotting the hard color of 
each source against its photon flux. These quantities are 
measures of the physical quantities of interest, the intrin- 
sic spectral shape and the luminosity of a system. We 
also indicate in Figure 13 the expected hardness ratios and 
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Fig. 13. — Summary of the spectral and variability properties of the X-ray sources detected toward the Galactic center. For all of the 
sources in the sample, we plot the hard color against the photon flux. The hard color is defined as HR = (h — s)/{h + s), where h is the 
number of counts in a hard band from 4.7—8.0 keV, and s is the number of counts in a softer band from 3.3-4.7 keV. We were able to derive 
hard colors for 1848 out of 2357 sources. Sources that are only detected in the 4.7-8.0 keV band are given values of 1. Sources that are only 
detected in the 3.3—4.7 keV band are given values of —1. Sources detected in neither band are given hard colors of —1.1. Sources located at or 
beyond the Galactic center are indicated in blue, while foreground sources are indicated in red. We use solid lines to illustrate the hard colors 
and photon fluxes expected for thermal plasma with a range of temperatures and luminosities (bolometric for d = 8 kpc), and the dotted 
lines for power laws with a range of photon indices and luminosities (0.5-8 keV). We also indicate which sources exhibit iron lines (stars) 
and short-term variability (diamonds) or long-term variability (circles). Note that sources with short-term variability were not searched for 
long-term variability. 



photon fluxes for sources over a range of luminosities and 
with either (1) thermal plasmas over a range of tempera- 
tures kT (solid lines), or power laws over a range of photon 
indices F (dotted lines). In all cases, we have assumed the 
sources are absorbed by the median column density from 
the spectral models, 6 x 10^^ cm~^ of gas and dust. Ap- 
proximately 75% of the 2000 Galactic center X-ray sources 
are detected with 90% confidence in both the 3.3-4.7 and 
4.7-8.0 keV bands, and therefore have hard colors in Fig- 
ure 13. 

We also indicate in Figure 13 which sources exhibit line 
emission from Fe. Nearly all of the sources brighter than 
4 X 10~^ ph cm~^ s~^ and harder than HR = exhibit 
line emission from He-like Fe. This is not surprising, given 
the prominence of line emission in the average spectra. Fe 
emission is detected less often in fainter sources, but this 
is probably due to lower signal-to-noise. 

Finally, we indicate which sources exhibit variability. 
The sources that are identified as variable tend to be 
brighter, because the signal-to-noise is better. Soft sources 
are most likely to exhibit short-term variability. 

In the following sections, we use these properties to 
guide our discussion of the natures of the sources. The 
luminosities of the Galactic center sources are consistent 
with those of young stellar objects (YSOs), interacting bi- 



naries (RS CVns), Wolf-Rayet (WR) and early O stars, 
cataclysmic variables, quiescent black hole and neutron 
star X-ray binaries, and possibly the ejecta of recent su- 
pernova that are interacting with molecular clouds. We 
consider each in turn. 

3.1. Sources with Active Stellar Coronae 

Many stars produce X-rays in their magnetic coronae. 
In particular, K and M dwarfs are so numerous that they 
contribute significantly to heating the ISM (e.g., Schlick- 
eiser 2002). However, individually their X-ray emission is 
faint, with Lx < 10^^ erg s~^, and cool, with kT < 1 keV 
(e.g., Krishnamurthi et al. 2001). Although, most of the 
foreground sources are probably low-mass main sequence 
stars (P. Zhao, in preparation), few of the Galactic center 
sources should be. YSOs and RS CVns are significantly 
brighter, with Lx ~ 10^^ to lO'^^ erg s^^(e.g., Feigelson 
et al. 2002; Dempsey et al. 1993a), and are therefore more 
likely to be seen at the Galactic center. 

3.1.1. Young Stellar Objects 

The number of YSOs at the Galactic center will depend 
upon whether low-mass stars have formed there recently. 
For instance, if star formation proceeds at the Galactic 
center in a similar manner as it has in the Orion nebula, 
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then the two-dozen massive, emission-hne stars in the cen- 
tral parsec of the Galaxy could conceivably be accompa- 
nied by tens of thousands of low-mass YSOs (e.g., Feigel- 
son et al. 2002). However, the strong tidal forces, mil- 
liGauss magnetic fields, and turbulent molecular clouds 
near the Galactic center may prevent low-mass stars from 
forming there (Morris 1993). 

YSOs have luminosities between 10^^ — 10^^ '' erg s~^, 
and spectra that can be described by thermal plasma emis- 
sion with kT =1 — 10 kcV (e.g., Pricbisch & Zinncckcr 
2002; Kohno, Koyama, & Hamaguchi 2002; Feigelsonet al. 
2002). Therefore, any YSOs that are located near the 
Galactic center should be found in the bottom left of Fig- 
ure 13, with HR < and fluxes < 2 x 10^'' ph cm^^ 
s""'^. However, the detection threshold for sources with 
HR < is approximately 10^^'^ erg s~^, and only ~ 0.4% 
of YSOs are brighter than this limit (Feigelson et al. 2002). 
YSOs also commonly exhibit flares lasting several hours, 
but fewer than 0.1% exhibit flares brighter than 10"^^ erg 
s~^ (Grosso et al. 2004; Feigelson 2004). In contrast, the 
faintest genuine flare from a Galactic center source has a 
peak luminosity of 5 x 10^^ erg s^^, whereas only three 
sources with short-term variability have peak luminosities 
below 10^^ erg s~^. Therefore, we believe that even the 
flaring sources are unlikely to be YSOs, and that any pop- 
ulation of YSOs remain largely undetected at the Galactic 
center. 

3.1.2. Interacting Binaries 

RS CVn systems are among the most numerous hard 
X-ray sources with Lx > 10^^ erg s~^, with a local space 
density of 5 x 10~^ pc~^ (Favata et al. 1995). Using 
the models of Launhardt et al. (2002) to scale the local 
number density to the stellar density at the Galactic cen- 
ter (Section 3), we estimate that the total number of RS 
CVns within 20 pc of the Galactic center is « 1.5 x 10^, 
while the number within a cylinder of 20 pc radius extend- 
ing the length of the nuclear bulge (440 pc) is 7 x 10*. 

However, RS CVns would be difficult to detect near the 
Galactic center. They typically have soft spectra, with 
fcT w 0.1 - 2 keV, and luminosities of Lx = 10^^ - 10^^ erg 
s~^(e.g., Dempsey et al. 1993b; Singh, Drake, & White 
1996). Therefore, RS CVns would have HR < -0.3 and 
photon flux < 3 x 10^^ ph cm^^ s^^ in Figure 13. This 
portion of the flgure is sparsely populated. Moreover, we 
are only sensitive to sources with kT < 2 keV if they are 
more luminous than 5 x 10^^ erg s"-'^, whereas only ^ 2% 
of RS CVns are this luminous (Dempsey et al. 1993a). 
Finally, although RS CVns do exhibit flares lasting sev- 
eral hours with amplitudes of up to a factor of ten, they 
are seldom more luminous than « 10^^'^ erg s~^ (e.g., 
Tsuru et al. 1989; Franciosini, Pallavicini, & Tagliaferri 
2001, 2003). Therefore, these flares would not be observ- 
able in our Galactic center data. The difflculty of detecting 
RS CVns and the lack of good candidate objects indicates 
we have probably identified only a tiny fraction of the RS 
CVns at the Galactic center. 

3.2. Winds from Massive Stars 

There is currently significant debate about the origin of 
the X-ray emission from WR and early O stars (e.g., Wal- 
dron & Cassinelli 2001), but it is generally thought that 



the X-rays are produced through shocks in their winds 
(see, e.g., Chlebowski & Garmany 1991). They have lu- 
minosities of up to « 1033-5 gj.g g-1 isolation, and 
» 10^^ erg s~^ when two such stars are in a colliding- 
wind binary. Their spectra can usually be modeled as 
thermal plasma with kT = 0.1 — 6 keV (e.g.. Pollock 1987; 
Portegies-Zwart et al. 2002). These systems would lie in 
the portion of the color-intensity diagram with HR < in 
Figure 13. 

The number of these systems present near the Galactic 
center is unknown, because it is determined by the uncer- 
tain star formation history (see, e.g., Morris 1993). Figer 
(1995) and Cetera et al. (1999) have identified several mas- 
sive, emission-line stars associated with HII regions near 
the Galactic center, but none of these has counterparts in 
our X-ray catalog. The wide-area search that Figer (1995) 
conducted failed to turn up additional candidates. Still, 
the unique conditions at the Galactic center make it im- 
portant to understand the number of massive stars there, 
so we suggest that the relatively soft X-ray sources would 
serve as good targets for future searches for massive stars. 

3.3. Millisecond Pulsars 

Isolated millisecond pulsars typically produce X-ray 
emission from particles accelerated as they spin down, with 
Lx = 10^^ - 10^1 erg s"^ (Possenti et al. 2002). At these 
luminosities, millisecond pulsars would be undetectable at 
the Galactic center. However, Cheng et al. (2004) have 
predicted that the wind from a millisecond pulsar could 
produce Lx = 10"^^ — 10"^^ erg s~^ by interacting with 
dense regions of the ISM (n > 100 cm~^). They suggest 
that ^ 100 millisecond pulsars could be present in our 
field, although the number of detectable systems would de- 
pend upon the volume of dense gas at the Galactic center. 
Muno et al. (2004) have demonstrated that a large fraction 
of the inner 20 pc of the Galaxy is filled with hot (T ~ 10^ 
K), low density (n w 0.1 cm^"*), X-ray emitting plasma, so 
only a small fraction of isolated millisecond pulsars may be 
detectable. Moreover, their spectra should be power laws 
with r = 1.5-2.5, which corresponds to -0.2 < HR < 0.0 
in Figure 13. This places millisecond pulsars on the same 
portion of the hardness-intensity diagram as CVs and RS 
CVns. As discussed in Cheng et al. (2004), identifying 
candidate systems among the point sources would be dif- 
ficult, but millisecond pulsars could account for extended 
features seen in the field (Morris et al. 2003). 

3.4. Accreting Sources 

3.4.1. Low-Mass X-ray Binaries 

Neutron stars and black holes accreting from low-mass 
companions that over-fill their Roche lobes are typically 
identified in outburst with Lx > 10^® erg s~^, although 
the majority of their time is spent in quiescence with 
Lx < 10'^^ erg s^-'^. LMXBs have been observed exten- 
sively in quiescence. The spectra of quiescent neutron star 
systems have been described with a fcT w 0.3 keV black 
body producing Lx ^ 10'^^ erg s^^, plus a F « 1 — 2 
power-law tail that contributes Lx ^ 10'^^ erg s^^ (Asai 
et al. 1998; Kong et al. 2002a). The black hole systems 
have Lx < 10^^ erg s~^ and exhibit F « 1 — 2 power- 
law spectra (Rutledge et al. 2001; Wijnands et al. 2002; 
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Campana et al. 2002a). The thermal emission from a neu- 
tron star would be unobservable behind 6 x 10^^ cm~^ of 
absorption, so is of little relevance to the current observa- 
tions. The power-law components of both the neutron star 
and black hole systems would produce —0.1 < HR < 0.2. 

However, LMXBs are rare — theoretical models pre- 
dict that ~ 10^ should currently be in quiescence in the 
entire Galaxy, whereas only only ~ 100 LMXBs, or 1%, 
have been identified (compare Iben, Tutukov, & Fedoroval 
1997; Belczynski & Taam 2004). Thus, if LMXBs form at 
the Galactic center in a similar manner as in the disk, our 
observation should encompass ~ 20 of them (Belczynski & 
Taam 2004). Transient outbursts from three LMXBs al- 
ready have been identified within 10' of the Galactic center 
(Eyles et al. 1975; Pavlinsky, Grebenev, & Sunyaev 1994; 
Maeda et al. 1996). If these truly represent ~ 20% of 
the total number there, then it would appear that LMXBs 
near the Galactic center are considerably more active than 
those in the Galactic disk. Alternatively, LMXBs could be 
concentrated near the Galactic center through dynamical 
settling (e.g., Morris 1993; Portegies-Zwart, McMillan, & 
Gerhardet al. 2003). In order to better constrain the num- 
bers of LMXBs within the nuclear bulge, it is important 
to continue to monitor this region in order to search for 
transient outbursts from additional systems. 

Prior to the Roche-lobe overflow phase, accretion onto 
the compact objects should also proceed at low rates 
from the winds of the low-mass companions (Bleach 2002; 
Willems & Kolb 2003; Belczynski & Taam 2004). These 
pre-LMXBs should have Lx = 10^^ - 10^^ erg s~S and 
would probably resemble Roche-lobe overflow systems in 
quiescence. Up to 10^ systems could be present in the 
Galaxy, and ^ 20 — 100 in our image of Sgr A* (Willems 
& Kolb 2003; Belczynski & Taam 2004). 

3.4.2. High-Mass X-ray Binaries 

Neutron stars and black holes accreting from the winds 
of massive companions should be about as common as 
LMXBs, because although the massive companions have 
much shorter lifetimes, accretion can occur when the sep- 
arations between the binary components are much larger 
(~ 1 AU compared to ~ Rq; see Pfahl, Podsiadlowski, & 
Rappaport 2002). They could be particularly abundant 
near the Galactic center, because it appears that 10% of 
Galactic star formation is currently occurring within the 
nuclear bulge (Launhardt et al. 2002). Our observations 
encompass ~ 5% of the nuclear bulge, so it would not be 
unreasonable to assume that, of the ~ 10"^ HMXBs in the 
Galaxy, 0.1 • 0.05 • 10^ ~ 50 could be present in the field 
around Sgr A* (sec Pfahl et al. 2002). 

In both outburst and quiescence, black hole HMXBs 
generally resemble LMXBs, because their X-ray emission 
is produced entirely in the accretion flow. Neutron star 
HMXBs, on the other hand, usually look much different 
from their LMXB counterparts, because the neutron stars 
in the young, high-mass systems tend to be more highly 
magnetized {B > 10^^ Gauss). Neutron star HMXBs in 
outburst produce X-rays from shocks that form in the 
magnetically-channeled column of accreted material. At 
the location of the shocks, the accretion flow is optically- 
thick, so the resulting spectra are flat, and can be de- 
scribed with a r < 1 power law between 2-8 keV (e.g. 



Campana et al. 2001). Therefore, neutron star HMXBs 
should have HR > 0.1 in Figure 13. HMXBs also some- 
times exhibit line emission at 6.4 keV from fluorescent neu- 
tral material in the companion's wind, as well as weaker 
emission from He-like Fe at 6.7 keV that is produced by 
photo-ionized plasma in the wind. Although large equiv- 
alent widths have been reported from low-resolution mea- 
surements with gas proportional counters (e.g., Apparao 
et al. 1994), the few measurements of these lines with 
CCD resolution spectra indicate that they have equiva- 
lent widths < 100 eV (e.g., Nagase et al. 1994; Shrader 
et al. 1999). 

Since the strong magnetic fields around the neutron 
stars in HMXBs channel the accreted material onto the 
star's polar caps, the surest way to identify neutron star 
HMXBs is through periodic modulations in their X-ray 
emission. We have found that seven hard Galactic center 
sources in our field exhibit periodic variability (Muno et al. 
2003c). However, the periods are all > 300 s, which makes 
it impossible to rule out that they are accreting white 
dwarfs. On the other hand, modulations with shorter peri- 
ods would be rendered undetectable by Doppler shifts from 
orbital motion, so the lengths of the periods observed are 
not necessarily a strong constraint on the entire population 
of sources at the Galactic center. 

Other variability is also seen from HMXBs. Short-term 
(several ks) flares are seen infrequently and are ascribed to 
instabilities in the accretion flow (e.g., Augello et al. 2003; 
Moon, Eikenberry, & Wasserman 2003). Long-term vari- 
ations are more common and are often caused either by 
changes in the density of the wind at the location of com- 
pact objects that have eccentric orbits around the donor 
star, or by instabilities in the excretion disks around the 
Be stars that are the mass donors in half of the known 
HMXBs (e.g. Apparao et al. 1994). 

The above considerations suggest that faint, neutron 
star HMXBs can account for some fraction of the hard 
Galactic center point sources. The main problem with 
this hypothesis is that few HMXBs have been observed 
at Lx < 10^^ erg s~^, and the ones that have can be 
described with much softer F ^ 2 spectra (e.g., Campana 
et al. 2002a). Nonetheless, the physics of X-ray production 
at low accretion rates is uncertain, so we cannot be certain 
of what X-ray properties to expect from faint HMXBs. 

3.4.3. Cataclysmic Variables 

Cataclysmic variables (CVs) are the most numerous 
accretion-powered X-ray sources. Their local space den- 
sity is ~ 3 X 10~^ pc"^ (Schwope et al. 2002), so that if 
we scale their number to the stellar density at the Galactic 
center, within 20 pc of Sgr A* we would expect ^ 9 x 10^ 
CVs, and within a cylinder centered on the Galactic center 
that is 20 pc in radius and 440 pc deep we would expect 
~ 4 X 10*. About 50% of CVs are luminous enough to be 
observed from the Galactic center (Verbunt et al. 1997), 
so they could account for the majority of the X-ray sources 
detected there. 

Systems with non-magnetized white dwarfs, which com- 
prise 80% of CVs, have luminosities between 10^^'^ — 
20^2 gj-g g-i^ and spectra that can be described with 
kT = 1 — 25 keV plasma from an accretion shock (e.g., 
Eracleous, Halpern, & Patterson 1991; Mukai & Shiokawa 
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1993; Verbunt et al. 1997). Thus, these systems should 
have hard colors HR < in Figure 13, and would be lo- 
cated in a similar portion of the color-intensity diagram as 
RS CVns and YSOs. 

CVs containing magnetized white dwarfs, which are 
referred to as polars and intermediate polars depending 
on whether or not the rotational period of the white 
dwarf is synchronized to its orbital period, comprise about 
20% of all CVs (e.g., Warner 1995, see also the CVcat 
database^''). Polars have similar spectra and luminosities 
as un-magnetized CVs, with the addition of a fcT ~ 50 eV 
"soft excess" that is attributed to "blobs" of accreted ma- 
terial that penetrate deeply into the photosphere (e.g., 
Ramsay et al. 1994; Verbunt et al. 1997; Ezuka & Ishida 
1999). The soft component woidd be unobservable above 
2 keV, so polars should also have HR < 0.1 in Figure 13. 
Polars also commonly exhibit variations in their average 
luminosity on time scales of years: « 50% of the polars 
surveyed by Ramsay et al. (2004b) changed in intensity 
by factors of > 4 between observations taken with ROSAT 
(1990-1999) and XMM-Newton (2000-present). Such vari- 
ations would be detectable from most of the Galactic 
center sources (Figure 11). Therefore, in the two years 
spanned by the Chandra observations of the Galactic cen- 
ter, we would expect ~ 15% of the polars to exhibit long- 
term variations. Since only 2% of the sources located at 
or beyond the Galactic center are variable, at most 20% 
could be polars. 

The intermediate polars are typically more luminous 
than other CVs, with Lx = lO^^ - lO^^-^ erg s'S and 
represent about 5% of the total population (see CVcat; 
Kube et al. 2003). This is thought to be related to the 
fact that they tend to have longer orbital periods (> 2 h), 
which could result in a higher mass transfer rate; however, 
the high M could also be a selection effect, because if a CV 
is bright, it is easier to detect modulations in the X-ray 
and optical emission at the rotational and orbital periods 
(Warner 1995). Intermediate polars also typically have 
much harder spectra than other CVs: when approximated 
as a power law, the optically thin thermal plasma usually 
seen from CVs should have F w 1.5, whereas the spectra of 
intermediate polars usually have F « 0. This is probably 
a result of the geometry of the accretion flow, because, as 
in other CVs, prominent line emission from He-like and 
H-like Fe indicates that the X-rays are produced either by 
plasma with kT « 1 — 20 keV or by a plasma photo- ionized 
by continuum X-rays that are not observed directly (e.g., 
Ezuka fc Ishida 1999; Mukai et al. 2003). In either case, 
the X-ray emitting regions would have to be partially ab- 
sorbed by material in the accretion flows, which removes 
low-energy photons from the spectra, thus making them 
flatter. Intermediate polars should have HR > 0.1 in Fig- 
ure 13, which makes them the best candidates among CVs 
for the hard Galactic center sources. 

The detailed spectral properties of intermediate polars 
are broadly consistent with the average spectra of the point 
sources in Figure 7. Weak emission at 6.4 keV is observed 
from these systems, and is attributed to X-rays that reflect 
off of the white dwarf's surface (e.g., Mukai & Shiokawa 
1993; Ezuka & Ishida 1999). Moreover, when the spectra 
of intermediate polars are modeled as emission from ther- 

http://minerva.uni-sw.gwdg.de/cvcat/tpp3.pl; Kube et al. (2003) 



mal plasma, the derived Fe abundances are often near or 
below the solar values (e.g.. Done & Osborne 1997; Fuji- 
moto & Ishida 1997; Ishida et al. 1997). This is similar to 
what we infer for the point sources in Table 5. 

Finally, the general lack of variability in the X-ray emis- 
sion from the Galactic center sources (aside from peri- 
odic modulations) is also consistent with the stable emis- 
sion usually seen from intermediate polars. On long time 
scales, the optical luminosity of intermediate polars usu- 
ally remains constant for many decades (e.g., Garnavich & 
Szkody 1988); because the optical and X-ray flux are cor- 
related in polars, we would expect the X-ray emission from 
intermediate polars also remains constant. Flares lasting 
several hours, presumably from accretion events, are some- 
times observed from magnetic CVs, but appear to be rare 
and most prominent in the soft X-ray band (< 2 keV; e.g., 
Patterson & Szkody 1993; Choi, Dotani, & Argawal 1999; 
Still & Mukai 2001). The predominant short time scale 
variability in intermediate polars is due to modulations 
of the emitting regions as the white dwarfs rotate (e.g., 
Norton & Watson 1989; Schwope et al. 2002; Ramsay & 
Cropper 2003). We have detected periodic modulations 
from seven of the brightest 285 Galactic center sources 
(Muno et al. 2003c). Since we were only sensitive to high- 
amplitude modulations, it is likely that many sources with 
low-amplitude modulations went undetected. Therefore, 
although the faintness of the Galactic center X-ray sources 
is probably the main cause of the lack of observed short- 
term variability, it is also plausible that the sources are 
intrinsically steady X-ray emitters like intermediate po- 
lars. 

Since the properties of the Galactic center sources 
change little as a function of their luminosity between lO '^ 
and 10'^^ erg s~^ (Figures 8 and 13), we believe that the 
majority of the Galactic center sources are intermediate 
polars. Intermediate polars comprise 5% of all known CVs 
(Kube et al. 2003), so given that there could be 4x lO'^ CVs 
within a pencil-beam centered on the Galactic center that 
is 20 pc in radius and 440 pc deep, they could reasonably 
account for the 1000 X-ray sources with HR > 0. 

3.5. Supernova Ejecta 

Bykov (2002, 2003) has suggested that the point sources 
in the Galactic center may not be stellar, but could be 
iron-rich fragments of supernova explosions that are inter- 
acting with molecular clouds. On order 10^ X-ray emitting 
knots could plausibly be produced by just 3 supernova oc- 
curring within the last 1000 y within 20 pc of the Galac- 
tic center; already, Sgr A East (Maeda et al. 2002) and 
the radio wisp 'E' (Ho et al. 1985) are thought to be rem- 
nants of recent supernova. The observational properties of 
the point sources can be reproduced by choosing several 
parameters in the ejecta model (Bykov 2003): the slope 
log AT — logS distribution of the knots (a » 1.7) is de- 
termined by their sizes and velocities, the slopes of their 
continuum X-ray emission (0 < F < 1.5) is set by the 
amplitudes of magneto-hydrodynamic turbulence in the 
shocks they produce, and the equivalent widths of the Fe 
emission (up to 1 keV) by their iron abundances. Future 
observations of known supernova remnants will better con- 
strain the properties of the X-ray emitting knots, which in 
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turn could make it possible to distinguish such knots from 
the stellar sources in the field. 

3.6. Unusual Sources 

A handful of the Galactic center sources resemble un- 
usual objects that have been found through shallower 
ASCA, BeppoSAX, XMM-Newton, and INTEGRAL sur- 
veys of the Galactic plane. These sources are important, 
because they could represent stellar remnants that are in 
short-lived states of accretion. We list the properties of 14 
unusual sources from other surveys in Table 8. The first 
three are polars that were identified with ASCA as having 
unusually strong emission lines from He-like Fe (equivalent 
widths > 1 keV) ; the fourth XMM-Newton source has sim- 
ilarly strong Fe emission at 6.7 keV, but its nature is un- 
certain. We find that 6 out of 183 Galactic center sources 
searched for Fe emission have 6.7 keV lines with equivalent 
widths greater than 1 keV, which is similar to the fraction 
of such sources identified in the ASCA Galactic plane sur- 
vey. The next four are highly- absorbed (A^h > 10^^ cm~^) 
sources identified with INTEGRAL and XMM-Newton, 
one of which has strong low- ionization Fe emission with an 
equivalent width > 1 keV. Wc find that 30% of the Galac- 
tic center sources have similarly high absorption, and two 
systems exhibit 6.4 keV Fe lines with equivalent widths > 1 
keV (CXOGC J174613. 7-290662 and GXOCG J174617.2~ 
285449 in Table 3). The final five are hard X-ray sources 
with slow (> ICQ s), high-amplitude periodic modulations 
in their X-ray emission. We find seven hard sources near 
the Galactic center (and one foreground source) with sim- 
ilar periodic X-ray modulations (Muno et al. 2003c). 

These sources would have been difficult to identify with 
the soft X-ray detectors on ROSAT (0.1-2.4 keV), which 
was the last observatory that systematically surveyed the 
sky for faint X-ray sources. Our study of the Galactic cen- 
ter suggests that they account for a few percent of all faint 
X-ray sources. 

4. CONCLUSIONS 

We have established that, on average, the X-ray sources 
detected in 626 ks of Chandra ACIS-I observations of the 
field around Sgr A* have hard, F < 1 spectra with promi- 
nent emission from He-like Fe at 6.7 keV (Figure 7 and 



Table 4). They also generally do not vary by more than 
factors of a few on time scales of hours or months. The 
best candidates for these hard X-ray sources are interme- 
diate polars, which represent the most luminous and spec- 
trally hardest 5% of all CVs. Therefore, the Galactic cen- 
ter X-ray sources are likely to be only a sub-sample of a 
population of ^ 10* CVs located near the Galactic center. 

Although a single population of sources may dominate 
the image, there are certainly many classes of objects 
present in smaller numbers in the field. Determining the 
numbers of rare objects is particularly important. For in- 
stance, the numbers of massive Wolf-Rayet and O stars 
and faint neutron star high-mass X-ray binaries can con- 
strain the recent rate of massive star formation near the 
Galactic center, while the numbers of LMXBS provide di- 
rect tests of the validity of unusual pathways for binary 
stellar evolution. For this reason, we are carrying out 
deep infrared observations of the Galactic center to iden- 
tify counterparts to the X-ray sources. These observa- 
tions will be useful for distinguishing CVs from, for ex- 
ample, HMXBs and WR/0 stars. At at a distance of 
8 kpc and with an extinction of Ak ~ 5 (Tan fc Draine 
2003), CVs should have K magnitudes of 22 25, and there- 
fore would be among the faintest detectable sources at the 
Galactic center (Warner 1995; Hoard et al. 2002). In con- 
trast, HMXBs and WR/0 stars should have K magnitudes 
brighter than 15 (Zombeck 1990; Wegner 1994) and will 
be very easy to detect. Therefore, the prospects for identi- 
fying the natures of the Galactic center X-ray sources are 
promising. 
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Table 1 

Observations of the Inner 20 pc of the Galaxy 









Aim Point 




Start Time 


Sequence 


Exposure 


RA 


DEC 


Roll 


(UT) 




(s) 


(degrees 


J2000) 


(degrees) 


1999 Sep 21 02:43:00 


0242 


40,872 


266.41382 


-29.0130 


268 


2000 Oct 26 18:15:11 


1561 


35,705 


266.41344 


-29.0128 


265 


2001 Jul 14 01:51:10 


1561 


13,504 


266.41344 


-29.0128 


265 


2002 Fob 19 14:27:32 


2951 


12,370 


266.41867 


-29.0033 


91 


2002 Mar 23 12:25:04 


2952 


11,859 


266.41897 


-29.0034 


88 


2002 Apr 19 10:39:01 


2953 


11,632 


266.41923 


-29.0034 


85 


2002 May 07 09:25:07 


2954 


12,455 


266.41938 


-29.0037 


82 


2002 May 22 22:59:15 


2943 


34,651 


266.41991 


-29.0041 


76 


2002 May 24 11:50:13 


3663 


37,959 


266.41993 


-29.0041 


76 


2002 May 25 15:16:03 


3392 


166,690 


266.41992 


-29.0041 


76 


2002 May 28 05:34:44 


3393 


158,026 


266.41992 


-29.0041 


76 


2002 Jun 03 01:24:37 


3665 


89,928 


266.41992 


-29.0041 


76 



Table 2 

Spectra of Individual Point Sources 









Power Law 






Plasma 


mekal) 




Fe Lines 


Name 


Not Nh 


r 


Fx tiFx 






kT Fx 


uFx 




6.4 keV 6.7 keV Flags 


CXOGC J 


Counts 10^^ cm-^ 




(10"" erg cm"^ 8"^) 




10^^ 


(keV) (10 


erg cm~^ s~^) 





174521.9- 


-290519 


174521.9- 


-290616 


174522.3- 


-290322 


174522.9- 


-285718 


174522.9- 


-290706 


174523.1- 


-290205 


174523.2- 


-290116 


174523.3- 


-290637 


174523.4- 


-290248 


174523.8- 


-290514 


174523.8- 


-290652 


174524.0- 


-285947 


174524.1- 


-285845 


174524.7- 


-290038 


174525.1- 


-285703 



499 
366 
97 
139 
158 

123 
93 

127 
80 
92 

94 
82 
224 
121 
152 



13 



10+_l° 

< 8 
3 

16 

6 

< 12 

< 10 

(,4-46 

< 26 

< 4 

3 

< 0.5 



0.1 
0.3 



+ 0.9 
-0.5 
4-1.6 
-0.9 

-1 0+^-^ 
-0.3(e) 
3.3+? " 



1.4 



-2.7 

+ 3.0 
2.1 

1.3(c) 

+ 1.9 
1.2 
+ 1.2 
1.1 



-1.6 
-1.9 
> -0.1 



-2.1 
-0.9 
2.9 



+ 5.S 
0.9 
+ 1.8 
-1.0 
+ 0.6 
-0.4 



-0.6(e) 

1 7+0.4 
^- ' -0.3 



4.0 
3.3 
0.8 
0.9 
0.8 

0.8 
0.5 
1.1 
0.7 
0.6 

0.9 
0.6 
0.2 
0.9 
0.2 



6.5 
5.9 
0.9 
1.1 
5.2 

2.0 
0.9 
1.2 
0.7 
1.1 

0.9 
0.7 
0.2 
1.0 
0.2 



24/27 
14/20 
6/4 
15/7 
14/16 

9/5 
7/3 
9/13 
1/3 
9/7 

15/13 
2/2 
2/9 
11/4 
4/6 



18+ 
16l 



8tl 

25l- 
6 

10 
< 



+ 16.1 
1.9 



63.2(e) 
15.6(e) 

> 1.6 
79.9(e) 

3.3' 

> 1.9 

> 1.7 

> 1.3 
79.9(e) 

> 0.9 
79.7(e) 

2 r + i.o 

6.0(e) 

r 9 + 9.3 



3.6 
3.1 
0.6 
0.8 
0.9 

0.8 
0.5 
0.9 
0.6 
0.6 

0.7 
0.5 
0.2 
0.7 
0.2 



9.7 
10.6 
2.0 
1.5 
6.5 

3.4 
1.1 
3.9 
2.8 
1.1 

3.7 
0.8 
0.2 
1.7 
0.2 



27/27 
15/20 
7/4 
19/7 
10/16 

5/5 
6/3 
12/13 
3/3 
9/7 

14/13 
6/2 
5/9 
13/4 
5/6 




Note. — Spectral parameters are marked with (e) when the uncertainty calculation failed to converge. Fx refers to the observed 2—8 kcV flux, and uFx refers to the 
dc-absorbcd 2—8 keV flux. The "Flags" column indicates those sources for which the spectra may be siispect, because they exhibit variability on the long-term (1) or short- 
term (s), lie near the edge of a CCD (e), or may be blended with a nearby source (c). The table in the printed version is only a sample of the full table, which is 
available through the electronic edition. The spectra, response functions, effective area functions, background estimates, and event lists for each source are available from 
littp: //www . astro . psu . edu/users/niel/galcen-xray-data/galcen-xray-data . html. 



Table 3 

Iron Emission from Sources with > 160 Net Counts 



Is3 

to 



Source 


r 


Low-ionization Fc (6.4 keV) 
(10"'' ph cm"^ s^i) (eV) 






r 


Hc-likc Fo (6.7 keV) 
/fc EWf^ 
(10^^ ph cm"^ s-^) (eV) 


Ax' 




xV'^ 


174508.7-290324 


1+3.6 


< 0.2 


< 830 


0.9 


0.064 


11.8/ 11 


3.4lt;*7 


r, 0+0.5 
"J-^-0.3 


1581 


3.3 


0.000 


9.3/ 11 


174510.3-285435 


3.0+Jd 


< 0.5 


< 1425 


1.2 


0.100 


21.6/ 17 


3.oll;° 


2+"-^ 


732 


2.0 


0.007 


20.7/ 17 


174510.5-290645 




0.9l°„-J 


411 


7.7 


0.000 


42.4/ 39 


1-61?;^ 


1-Oio5 


474 


8.6 


0.000 


41.2/ 39 


174512.4-290604 


1 n-i-3 4 


< 0.5 


< 1447 


4.2 


0.024 


11.0/ 20 


n + 3 1 

2-2I2.0 


n 3+"-'^ 


1085 


5.4 


0.000 


10.2/ 20 


174517.3-290440 


n n + O 6 

-0-3Io.l 


< 0.4 


< 558 


5.3 


0.011 


21.6/ 23 


I-. 6 

-o-2lS;4 


n ^ +0 3 

"-^-0.2 


467 


6.7 


0.002 


20.0/ 23 


174519.8- 290114 

174520.9- 285818 


— 9 

4 2+1- = 

^•^-0.7 


0.3l°;l 


> lO'' 


6.1 


0.003 


36.5/ 15 
10.3/ 11 


1 4+-1-* 

— 1.4 

4.8ll;°3 


c; + -0.5 
3+0-1 


948 
> lO'' 


4.9 
5.9 


0.006 
0.003 


28.3/ 15 
10.6/ 11 


174525.5-290028 


^■■^-0.5 


< 0.6 


< 449 


4.8 


0.012 


49.9/ 28 




^■"-0.4 


818 


13.0 


0.000 


33.0/ 28 


174527.6-285258 






223 


6.2 


0.009 


48.2/ 40 


-. 6 
1 -210.8 


q+"-^ 


310 


10.0 


0.000 


42.9/ 40 


174527.8-290542 


^■"-1.5 


< 0.2 


< 394 


1.1 


0.152 


6.7/ 13 


2.ol°;5 


^ 1 +0 4 
O-4I0.4 


761 


4.2 


0.003 


5.1/ 13 


174529.0-290406 


„ (.+0.6 
* — 1 4 


< 0.5 


< 699 


2.9 


0.028 


13.7/ 11 


2-0li ,^ 




716 


4.7 


0.006 


11.1/ 11 


174529.6-285432 


1.5i- = 


2+"-^ 


238 


3.3 


0.007 


11.4/ 18 


1^-^-0.8 


< 0.4 


< 642 


1.8 


0.050 


12.5/ 18 


174531.1-290219 


-1-51-1.5 








14.0/ 5 


^•"-1.9 




2450 


4.8 


0.003 


3.4/ 5 


174532.3-290251 


^ ^-1-0 7 
^•"^-1.5 


< 0.5 


< 345 


2.1 


0.060 


24.9/ 14 


-1.3 




589 


7.5 


0.000 


14.5/ 14 


174532.4-290259 


n A' 

^•^-1.0 


« .4-0 2 
"•^-0.3 


434 


4.9 


0.001 


17.5/ 12 


n J 4-1 2 


o.5l°;^ 


542 


6.5 


0.002 


14.1/ 12 


174534.5-285523 




< 0.2 


< 280 


0.9 


0.228 


22.2/ 16 


1.31^4 


Q 7+0.3 
' -0.3 


943 


12.3 


0.000 


6.5/ 16 


174534.5-290201 
174534.9-290118 


7+°'* 

^- ' -0.9 


< 0.5 


230 

< 321 


8.4 
2.6 


0.001 
0.053 


43.7/ 39 
20.0/ 17 


4+°-^ 

"-^-0.6 

1 O+0.7 
1-^-1.4 


6+°-^ 


344 
439 


14.3 
7.5 


0.000 
0.001 


35.7/ 39 
13.7/ 17 


174535.6-290034 




< 0.2 


< 122 


1.0 


0.148 


26.2/ 24 


, .-4-0 9 
J^-^-1.1 


"•"-0.4 


513 


7.5 


0.000 


19.1/ 24 


174536.1-285638 


n A-i 3 4 










259.9/111 


n ->4-0 1 
3-6I0.2 


9 7+0.4 
'-0.4 


2173 


48.5 


0.000 


156.4/111 


174537.6-290144 


4+" '' 

^ — 5 


< 0.4 


< 410 


3.5 


0.018 


23.5/ 21 


4+0-6 
- —0 5 


2+"-^ 


249 


3.7 


0.007 


23.3/ 21 


174537.7-290002 


-0 0+"-* 


< 0.4 


< 815 


4.8 


0.020 


22.7/ 15 




o-5lS:2 


890 


6.7 


0.000 


18.9/ 15 


174537.9-290134 




< 0.3 


< 15803 


3.9 


0.010 


12.2/ 9 


^-"-0.9 


0.21°;^ 


16521 


3.9 


0.007 


12.3/ 9 


174540.1-290055 


9 9 + 0.8 


0.2l°.| 


130 


3.4 


0.000 


41.1/ 51 


2 l+O-''' 


< 0.2 


< 129 


1.0 


0.050 


43.1/ 51 


174540.5-285550 


k4-0 9 

0-5Io.3 


< 0.4 


< 249 


1.7 


0.092 


42.2/ 26 


rt-t-O 4 

1-0I5.9 


Q 7+0.3 
"■'-0.3 


448 


8.0 


0.002 


31.7/ 26 


174541.2-290210 


-1 0+^-° 


< 0.7 


< 279 


2.1 


0.055 


26.1/ 17 




8+°-^ 

"■°-0.4 


361 


6.4 


0.001 


19.1/ 17 


174541.5-285814 




3+°-^ 

U.a_Q 2 


174 


6.7 


0.003 


78.1/ 78 


0.9t.°i 


< 0.5 


< 271 


3.9 


0.020 


81.1/ 78 


174541.6-285952 


-0 1+°'^ 


4+°* 


365 


5.1 


0.003 


12.2/ 14 


r. 0+0.8 
"-^-0.8 


< 1.1 


< 1060 


2.3 


0.021 


15.7/ 14 


174541.8-290037 




< 0.7 


< 432 


4.6 


0.019 


41.5/ 27 


0.3+1-5 


^•"-0.4 


678 


11.1 


0.001 


30.0/ 27 


174542.0-285824 


0.31- 


< 0.4 


< 451 


1.6 


0.112 


8.8/ 6 


0-3lli 




584 


4.4 


0.002 


4.3/ 6 


174542.2-285732 


-0 2+°-'' 




229 


4.1 


0.007 


23.0/ 18 




< 0.5 


< 601 


4.9 


0.020 


21.8/ 18 


174543.3-285605 




< 0.2 


< 342 


1.7 


0.083 


9.6/ 12 




n 5+°-* 


782 


8.3 


0.000 


4.0/ 12 


174543.4-285841 


1.2l°;t 


< 0.5 


< 312 


3.0 


0.029 


26.8/ 21 


1 4+0.6 


' -0.4 


409 


7.7 


0.000 


20.2/ 21 


174543.7-285946 
174544.3-290156 


0.711;? 


1 O + 0.5 
^••^-0.5 

2+°-^ 


654 
355 


14.1 
6.0 


0.000 
0.001 


23.9/ 28 
1.2/ 7 


7+0.6 
"•'-0.9 

"•^-1.9 


< 1.3 


493 
< 1890 


7.2 
1.8 


0.001 
0.060 


34.9/ 28 
3.8/ 7 


174544.8-285953 


o.3l?;° 


< 0.3 


< 1151 


0.9 


0.200 


32.1/ 17 




9+1^ 


5241 


8.7 


0.002 


17.5/ 17 



o 

£0 



Table :> Coiitmufd 







Low-ionization Fe (6.4 keV) 






He-like Fe (6.7 keV) 




Source 


r 






r 










(10-'' ph cm-^ s-^) (eV) 






(10"'' ph cm-= s-^) (eV) 





174544.9-290027 

174546.1- 290057 

174546.2- 285906 
174546.9-285903 
174547.0-285333 

174547.2- 290000 
174548.9-285751 

174549.3- 285557 
174549.6-290457 
174550.9-285430 

174552.0-285312 

174554.4- 285816 

174555.6- 285600 
174558.9-290724 

174559.5- 290601 

174601.0- 285854 

174601.1- 285953 

174606.3- 285810 

174608.4- 290623 

174609.8- 290321 

174610.9- 285345 

174612.3- 285706 

174613.7- 290622 

174614.5- 285428 

174616.5- 285846 

174617.2- 285449 

174619.4- 290213 

174623.6- 285629 



0.5 



-0.3 
+ 0.7 



0.4 



-1-7-0.5 
,+o.s 

-0.6 



0.31 



0.1 
0.! 



+ 2.3 
0.8 
+0.4 
-0.4 
+0.4 



0.91 



0.7ir4 

□+0.5 
-0.4 

,+0.3 



0.3] 
1.7 



0.3 
+0.6 
-1.9 



1 r+0.8 
^•''-0.9 

1 1+0.6 

„ „ + 0.5 
^■^-1.3 

^■"-2.6 

^••'-i.o 

1 3+^-^ 
^••5-1.1 

i+''-^ 

1 O+0.5 
-^■•5-0.5 

1.0 ' " " 



+ 0.2 
0.3 

< 0.3 

< 0.3 

< 0.4 

< 0.4 



0.5^ 



- + 0.3 
-0.2 

5 + '''2 
^•^-0.3 

^- ' -0.3 

-+0.3 
-0.4 



0.5: 
1.1^ 



1+0.8 
-0.5 

< 0.5 
1 1+°* 

< 0.5 



-0-5 

0.3 



< 1.1 

1 1 + "" 
-^•■'-0.5 

< 11.6 
7+0.6 

< 0.2 
0.3 
0.5 
0.3 
0.2 



+0.2 
0.2 
+ 2.1 
0.5 
+ 0.3 
0.3 
+ 0.4 
-0.2 



+0.5 
0.6 
+ 1.9 
-1.0 



7+°''^ 

' -0.5 

< 1.7 



237 

< 485 

< 292 

< 558 

< 188 

561 
261 

364 
414 
426 

238 

420 

< 1595 
338 

< 552 

< 284 
749 
841 

455 

< 1560 

336 

< 450 
1594 
923 
344 
3021 

369 

< 582 



6.7 
2.6 
3.1 
3.2 
1.4 

6.8 
7.7 
10.4 
6.5 
6.5 

6.9 

9.2 
3.4 
13.5 
3.0 

3.3 
4.5 

14.6 
9.9 
4.1 

4.8 
1.2 
6.8 
3.7 

5.5 
1.7 

7.6 
4.4 



0.000 
0.050 
0.034 
0.014 
0.124 

0.001 
0.000 
0.000 
0.000 
0.000 

0.001 

0.000 
0.020 
0.000 
0.050 

0.030 
0.006 
0.000 

0.000 
0.037 

0.001 
0.180 
0.006 
0.004 
0.006 
0.000 

0.001 
0.020 



33.0/ 38 
16.4/ 11 
21.6/ 31 
11.0/ 11 

60.5/ 43 

25.4/ 18 
66.4/ 49 
44.9/ 37 
12.2/ 14 
18.0/ 23 

59.1/ 43 

39.0/ 26 
6.9/ 12 
96.4/ 66 
29.6/ 20 

50.6/ 31 
4.2/ 7 
25.9/ 29 

48.3/ 36 
26.4/ 15 

48.4/ 51 
27.0/ 23 
15.3/ 22 
62.7/ 51 
23.2/ 27 
35.2/ 26 

13.1/ 22 
76.8/ 51 



+ 0.2 
-0.3 

-1 2+°-* 

' -0.4 

^■°-2.9 
7+0.4 
-0.4 



0.4 



0.7^ 



1 0+^'* 

1 + °'^ 

-0 2+°-^ 

' -0.4 
' -0.8 

-0 2+''-* 

5 5+"-^ 
"'•°-3.6 

'■'■^-0.3 

1 4+^'° 
-0 4+^-^ 

1 + °'^ 
n q + 0.8 
^•■^-1.4 

1 6+"-^ 
1 2+1-1 

^-1.4 
-1+1.2 
•^-1.3 

-i-4ii:^5 
o.9;S:'5 



0.4^ 



0.1^ 
1. 



1+0.2 
-0.4 

4+"-^ 
2+°-^ 
0.7+11 



0.. 



0.5 
+ 0.4 
-0.4 



4+1-5 
^-0.4 

1.6+?? 



-0.4'' 



0.8t°J 

< 0.8 

< 0.6 

8+°'* 

9+" '' 

1 4+°'* 
^•*-0.5 

4+°* 

"•^-0.3 

1 8+°* 

^•°-0.8 

< 0.4 

6+°'^ 

1 4+i'i 

^■^-0.7 

< 0.7 

'-'•'^-0.2 

< 0.5 

< 18.6 

< 0.4 

< 0.2 

< 0.9 

1 2+1" 

1-^-0.7 



248 
674 
217 
823 
423 

327 
446 
413 

< 862 

< 509 

183 
417 
2640 
416 
508 

526 

< 601 
458 
510 
623 

< 346 
817 

< 2798 

< 26553 

< 385 

< 2647 

< 432 
440 



3.9 
6.5 
5.8 
4.4 
10.8 

3.0 
13.8 
9.7 
2.8 
1.9 

3.9 
8.4 
5.5 
13.7 
5.0 

12.2 
0.9 
6.1 
8.9 
5.6 

1.8 
5.8 
1.7 
1.0 
1.3 
1.0 



0.001 
0.001 
0.000 
0.004 
0.000 

0.008 
0.000 
0.000 
0.013 
0.030 

0.002 

0.000 
0.003 
0.000 
0.002 

0.000 
0.170 
0.004 
0.002 
0.004 

0.100 
0.000 
0.060 
0.360 
0.160 
0.130 



1.6 0.090 
5.8 0.002 



35.8/ 38 
9.6/ 11 

19.6/ 31 
9.6/ 11 

47.3/ 43 

33.2/ 18 
56.8/ 49 
46.1/ 37 
17.5/ 14 
22.7/ 23 

63.8/ 43 

40.8/ 26 
5.5/ 12 
96.0/ 66 
26.3/ 20 

34.8/ 31 
8.5/ 7 
40.2/ 29 

50.0/ 36 
23.1/ 15 

51.5/ 51 
21.5/ 23 
20.1/ 22 
66.2/ 51 
27.6/ 27 
36.1/ 26 

18.2/ 22 
74.6/ 51 



Note. — Ppe refers to the chance probability that the line emission could be produce by random variations in a featureless continuum. Only sources with Ppf. < 0.01 
for a line at either 6.4 or 6.7 keV are listed in the table. Note that the line emission for those sources with a best-fit power laws with F > 4 should be viewed extremely 
skeptically, because these steep spectra are background-dominated above 6 keV. 
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Table 4 



Line Emission from Galactic Center Point Sources 





C < 500 


C < 80 


80 < C < 500 


HR > 0.1 


—0.1 < HR < 0.1 


HR < —0.1 


Nh (10^^ cm"^) 






9 9 + 0.2 
^■^-0 1 


1 4+0.2 

-0 1 


9 O + 0-3 


4 1 + 0-4 
■ -0 1 


Afpc.H (10^^ cm-^) 


7 1 +0.2 


7 n+O.l 


7 ,; + 0.5 


o O + 0.7 
8.3_(, 2 


o 9 + 0.8 


o 7 + 0.9 
'-0.4 


/pc 

r 


0.95 


0.95 


0.95 


0.95 


0.95 


0.95 






1 04+°"' 


34+0 12 


J-.OJ--0.04 


o 0-, +0.25 
■-'■■Jl-0.12 


Nr (10~® ph cm~^ arcmin"^) 


5.01°:? 


^••'-0.5 


o 9 + 0.8 
^■^-0.1 


1 9+0.3 


O-'-0.4 


29.2+3*2^ 


Si XIII He-a (eV) 


181 ± 29 


138 ± 38 


248 ± 57 


242 ± 65 


97 ± 41 


253 ± i23 


Si XIV Ly-a (cV) 


53 ± 14 


32 ± 20 


81 ± 26 


88 ± 34 


< 37 


85 ± 41 


Si XIII Ho-/j (cV) 


57 ± 15 


64 ± 24 


50 ± 24 


78 ± 35 


56 ± 27 


67 ± 42 


S XV Hc-a (cV) 


117 ± 19 


45 ± 20 


184 ± 33 


131 ± 35 


69 ± 26 


259 ± 128 


S XVI Ly-a (cV) 


21 ± 8 


< 18 


33 ± 13 


67 ± 20 


10 ± 13 


< 23 


Ar XVII Ho-a (cV) 


17 ± 5 


< 14 


23 ± 8 


35 ± 11 


< 4 


34 ± 18 


Ca XIX He-a (eV) 


7 ± 4 


5 ± 5 


7 ± 5 


< 2 


17 ± 7 


24 ± 14 


Fe K-a (eV) 


137 ± 21 


96 ± 21 


180 ± 32 


134 ± 27 


128 ± 37 


226 ± 120 


Fe XXV Hc-a (cV) 


404 ± 59 


465 ± 90 


350 ± 61 


396 ± 76 


388 ± 108 


411 ± 213 


Fe XXVI Ly-a (oV) 


225 ± 34 


266 ± 53 


195 ± 36 


209 ± 41 


217 ± 62 


335 ± 180 


F (10~^^ erg cm~^ s"^ arcmin"^) 


0.45 


0.22 


0.21 


0.26 


0.15 


0.04 


uF (10""^^ erg cm~^ arcmin"^) 


0.45 


0.22 


0.21 


0.26 


0.15 


0.04 




490/464 


448/464 


523/464 


593/464 


526/464 


518/436 



Note. — The line strengths are given in units of equivalent width. The fluxes are listed with units of arcmin ^ to indicate the 
amount of flux that is produced by point sources per unit area on the sky. 



Table 5 

Two-kT Plasma Models of Galactic Center Point Sources 





C < 500 


C < 80 


80 < C < 500 


HR > 0.1 


-0.1 < HR < 0.1 


HR < -0.1 


JVhi (10^^ 


4 5+"-" 


1-411:1 


4 8+"--'' 


1 2+=i 


;^ ;^+0.5 


7g+0.9 


JVpcHi (10^^ cm-^) 


^■°-0.4 


o 9 + 3.0 
"^■^-0.8 


4 3+1 = 

^■■^-0.8 


o -1 +0.8 
"^■1-1.0 


4 5+2-0 


22.81°:^ 


/pc,l 


0.95 


0.95 


0.95 


0.95 


0.95 


0.95 


feTi (keV) 


n r;f:+0.04 


67+° °* 

U.O' -0.08 


„ po+O.lO 


o.69i°o:;? 


69+° i° 


f, co+0.23 
'J-"°-0.09 


EMi (10"" cm"'^ pc) 


2 n+=° 




7+^ = 

' -0.4 


o.ol°o:^o 


0+°° 


' -5.1 


Fi (10~i^ erg cm"^ s~i arcmin"^) 


0.03 


0.01 


0.02 


0.01 


0.01 


0.01 


uFi (lO"!"^ erg cm~^ arcmin"^) 


0.2 


0.1 


0.1 


0.0 


0.0 


1.7 


JVh2 (10== cm-2) 


15.3l°:^2 


'•°-6.4 


16.21J:* 


8-5l?:J 


4 9+°-* 


^- ' -8.4 


JVpc,H2 (10== cm-=) 




24+101 


1191- 


39lt, 


10+1 


9+16 


/pc,2 


0.77 


0.66 


0.74 


0.81 


0.81 


0.95 


feT2 (keV) 


7 3+0.4 


9 0+°-^ 


7 g + 0.3 


8 4+° = 


8 6+°-* 

^'"-0.3 




EM2 (10-* cm-^ pc) 


1 71+0.47 
'■■ ' ^-0.24 




'O— 0.22 


44+° °* 

^■^^—0.02 


14+°°i 
"■i*-o.oi 


0.02+_°°\ 


F2 (10-1^ erg cm-= arcmin-=) 


0.41 


0.20 


0.18 


0.24 


0.15 


0.03 


UF2 (IQ-i^ erg cm-= s~' arcmin-=) 


3.60 


1.71 


2.10 


1.01 


0.32 


0.06 


Zsi/Zsi,Q 


'^■l''-0.05 


0.82lS;« 


26+°"^ 


1 64+°-''** 


r^7+0-58 


46+° ='' 




30+° "^ 




47+° i° 

^■^' -0.09 


1 t;-|+0-87 


1 14 + 0-79 


n 89+0-52 




45 + °-"-'> 


o.97lS:"ol 


n ,q+o.05 


67+° °* 


„ o-i +0.05 
0-O1-0.05 


f) 79+0.29 
l'-'''-0.02 


Fe K-a (10"^ ph cm-= arcmin-^) 


3-8-o'.8 




2 4+1.3 


1 4+°-i 


6+°-^ 


1+°° 




598/465 


564/465 


495/465 


622/465 


546/465 


520/437 



Note. — The fluxes are listed with units of arcmin"'^ to indicate the amount of flux that is produced by point sources per unit area 
on the sky. The abundances of Ar and Ca were held to their solar values. Note that the soft component of the spectrum does not 
improve the fit for the hard sources, with HR > 0.1 and —0.1 < HR < 0.1; we have included it only for the sake of comparison with 
the other groups of sources. 



Properties of Galactic Center X-ray Sources 



Table 6 

Sources with Short-term Variability 



oourcs 


~Loc — 


ObsID 


Var 
Type 


Al 
(ks) 


(10"' ph cm-^ s~^) 


(10"^ ph cm"^ s"^) 


F / F . 

^ max / ^ mm 


174517.4-290650 




3392 






28+^ 


105+" 


3 


1 745 2 . 3—290 1 43 




2943 


step 




0+5 




25 


174520.6-290152 


f 


3392 


flare 


91.0 


28+S 


1160+.?? 


40 






3393 


flare 


34.1 


16+^ 


819+S* 


51 


174521.8-285912 


f 


3392 


flare 


24.2 




1891=^ 


44 


J. i '-t'.) Z. '.) , A. Z^O^J 1 UO 


f 


3665 


step 






288lip 


28 


174530.3-290341 


gc 


3392 


flare 


22.0 




86I18 


14 


174531.0-285605 


gc 


3393 


step 




+ 2 
2-2 


40 9 


18 


174533.4-285328 


f 


0242 


flare 


8.2 


nr'-i-ll 

2619 


217l^J 


8 


i / 4004.0— ^yu^oo 


gc 


3392 


step 




< 1 


in+4 


> 5 


174535.6-290133 


gc 


3665 


flare 


61.2 


14+8 




11 


174535.8-290159 


gc 


3393 


step 




< 3 


12I* 


> 2 


174535.9-290806 


gc 


3393 


step 






3011 


6 


174536.3-285545 


f 


3392 


flare 


3.8 


0+2 


1821S 


63 


174538.2-285602 


f 


3665 


flare 


28.4 


i3i- 


116+^^ 


9 


174538.3-290048 


gc 


3392 


step 






5611 


2 


174540.1-290804 


f 


3392 


flare 


87.2 


< 3 


276lf, 


> 76 






3393 


step 




241- 


62l| 


2 


174540.4-285831 


gc 


3392 


step 






231^ 


4 


174541.4-290348 


gc 


3392 


step 




< 2 


20I5 


> 8 


174541.5-290752 


f 


2943 


flare 


7.6 


< 11 


124lfo 


> 7 


174541.8-290319'' 


gc 


1561a 










174542.8-285352 


gc 


0242 


step 




< 3 


361-, 


> 6 


174542.9-285522 


f 


3665 


flare 


3.3 


< 4 


1901^1 


> 30 


i 1 4040.4— ZyUo4 / 


gc 


1561a 










174543.9-290456 


f 


3392 
3663 


flare 
flare 


11.1 
6.1 


8ll? 


13971J?? 
433I77 


17 
3 


174545.0-290336 


gc 


3663 


step 




< 8 


37l- 


> 3 


174546.8-290252" 


gc 


1561b 










174547.4-290817 


f 


3393 


step 




' -4 


54l-a 


8 


174548.0-290352 


f 


3665 


flare 


11.0 




84lf, 


22 


174548.4-290234 


gc 


2943 


flare 


9.9 


< 22 


84lfe 


> 2 


174548.4-290832 


gc 


3393 


step 




< 3 


21I- 


> 4 


174548.6-290522 


f 


3392 


flare 


35.3 


12tl 


226^ 


18 


174550.7-290434 


f 


3393 


flare 


15.0 


< 4 


46+21 
^"-17 


> 7 


174552.1-290422 


gc 


3393 


step 




c+3 
"-3 


2611 


5 


174552.2-290744 


gc 


3392 


flare 


0.1 


99+^* 


2123'' 


21 


174552.9-290358 
174556.9-285819» 


f 
f 


0242 
3393 
2952 


flare 
flare 


0.7 

7.3 


341- 
201- 


652lfe°4 
40411- 


19 
20 


174558.5-290451 


f 


3392 


flare 


12.6 


1411 


253ir, 


18 


174559.0-290418 


f 


3392 


step 




< 3 


24l^e 


> 5 


174605.2-290700 


gc 


3393 


flare 


24.8 


< 4 


26911^ 


> 47 


174612.4-290234 


f 


3393 


step 




< 2 


17l^ 


> 3 



^The variability from these sources was identified via the KS-test, but not by the Bayesian Blocks algorithm. 
"^This flare consisted of 5 photons received within 100 s. 
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Table 7 

SonU'ES A\MTII LONC-TERM \'.VmABILITY 



Source 


Loc. 


ObsID 
of Min. 


(10-'^ ph cm"^ 


ObsID 
of Max. 


p 

(10-'' ph cm-2 s-^) 


J'max/J'min 


174503.9-290051 
174507 0— 29n'^5fi 


gc 


3392 
1561b 


441 6 
< 9 


0242 
2952 


00+24 

''*'-20 
41+25 


2 
> 2 


174514.1-285426 


gc 


2953 


< 13 


3392 


30+5 


> 1 


174517.5-285646 


gc 


3663 


< 3 


3665 


14+^ 


> 2 


174519.8-290114 


gc 


3663 




0242 




3 


174520.5-285927 


gc 


3665 


< 2 


2953 


00 + 29 

22I16 


> 2 


174520.6-285712 


gc 


0242 


< 8 


2954 


150+3'' 


> 15 


174520.8-285304 


f 


3665 


< 3 


3393 


15+4 


> 3 


174521.7-285812 


gc 


3665 


< 2 


2952 




> 1 


174521.9-290616 


gc 


1561b 


< 19 


2953 




> 2 


174522.4-285707 


gc 


3665 


< 2 


2953 


01 +22 
21-1.3 


> 3 


174523. 1-290205 


gc 


1561a 


< 3 


2954 


i,+13 




174526.4-290148 


gc 


3392 


< 1 


2951 


8+1^ 


> 1 


174526.7-290220 


gc 


3665 


< 2 


2952 


18+fo 


> 4 


174527.4-285938 


gc 


0242 


< 4 


2943 




> 4 


174529.0-290406 


gc 


1561b 


< 8 


3663 


25^^ 


> 2 


174529.6—285432 


gc 


2952 


< 9 


3393 


40+" 


> 3 


174530.5-290323 


gc 


3392 


< 1 


3663 


stl 


> 5 


174531.3-285949 


gc 


3393 


< 1 


2953 




> 6 


174531.8-290000 


gc 


3392 


< 1 


2951 




> 


174532.7-290552 


f 


0242 


< 3 


2954 


108+27 


> 28 


1 74532 Q— 285823 


gc 


3393 


4+2 


2952 


20+^" 


5 


174533.0-285355 


gc 


0242 


< 5 


3393 


38+5 


> 6 


174534.2-290119 


gc 


3393 


< 3 


2951 


74+11 


> 20 


174535.5-290124 


gc 


0242 


< 3 


2951 


25511^ 


> 80 


174535.9-285825 


gc 


3392 


< 1 


1561b 


9^7 


> 1 


174536 1—285638 
174536.6-290109 


f 

gc 


2952 
3393 


65+^** 

00_2i 

< 2 


1561b 
1561a 


iyy_3g 
13+f 


3 
> 2 


174537.2-285459 


f 


1561a 


< 3 


2951 


59+2I 


> 11 


174537.5-290125 


gc 


3665 


< 2 


2951 


24t\l 


> 5 


174538.0-290022 


gc 


3665 


39l6 


0242 


290tll 


7 


174538.4-290044 


gc 


1561a 


< 2 


2954 


108+3° 


> 33 


174538.7-290134 


gc 


1561a 


< 5 


3392 


19+3 


> 2 


174539.1-290112 


gc 


1561a 


< 4 


3392 


10+3 


> 2 


174539.5-285454 


f 


3392 


3l2 


2951 


25t\l 


8 


174540.1-290055 


gc 


2954 


I6I12 


3663 


oc + 16 


5 


174540.6-290001 


gc 


0242 


< 9 


1561b 


T 1 + 3 3 

134I28 


> 11 


174540.8-290040 


gc 


0242 


< 3 


2953 




> 


174541.0-290014 


gc 


1561a 


< 5 


3392 


204+S 


> 37 


174541.5-285148 


f 


3393 


< 5 


2953 




> 3 


174541.5-285814 


f 


2951 


< 14 


3663 




> 11 


174541.7-285555 


gc 


3663 


< 4 


3392 


241^ 


> 4 


174542.2-285732 


gc 


3665 


2 + % 

— 2 


3392 


57+« 

— 6 


25 


174542.2-290132 


gc 


3393 


< 1 


2951 


14"?/' 


> 4 


174542.5-285722 


gc 


3665 


< 1 


1561b 


10+;.! 


> 2 


174543.4-285742 


gc 


0242 


< 2 


2953 


37^30 


> 7 


174543.4-285900 


gc 


0242 


< 3 


2951 


0^+43 
8OI32 


> 15 


1 74543 6—285629 


gc 


3665 


0+3 
^-2 


2952 


0'3+45 


9 


174543.9-290245 


gc 


0242 


< 3 


2952 


I6+I0 


> 2 


174544.2-290644 


gc 


3665 


< 4 


3392 


13+3 


> 2 


174545.2-285828 


f 


1561b 


68l^? 


2951 


2Q8ttl 


3 


1 74546. 1-285831 


gc 


3665 


< 5 


1561b 


no + 21 

391 16 


> 4 


174546.6-290356 


gc 


3393 


< 1 


2953 


1 1 +13 


> 3 


1 74547 R— 2Qni 45 


gc 


1561a 


^ 5 


3663 


21+^ 
2^-7 




174549.5-285815 


gc 


3392 


< 2 


1561a 


lll^ 


> 2 


174550.5-285239 


f 


1561a 


O0 + 26 


1561b 


112^^4 


4 


174550.9-285430 


gc 


3393 


29t^ 


2951 


88l|3 


3 


174552.0-290324 


gc 


0242 


< 6 


2953 


22+_\l 


> 1 


174552.5-285759 


gc 


3393 


< 1 


2943 


io+_l 


> 4 


174553.3-290444 


f 


2943 


^tl 


3393 


46+_l 


8 


174553.3-290632 


gc 


3665 


itt 


2954 


ss+_ll 


10 


174554.2-285729 


gc 


3665 


2tl 


3392 


io+_l 


4 
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Table 7 — Continued 



Source 


Loc. 


ObsID 
of Min. 


Fmin 

(10"'' ph cm~^ s~^) 


ObsID 

of Max. 


Fmax 

(10"'' ph cm~^ s~^) 


.^max / Fr 


174558.4-290120 


f 


3665 


< 3 


1561a 




> 7 


174558.9-290724 
174601.0-285854 


gc 
gc 


3665 
2952 




0242 
1561b 


46i- 


2 
4 


174601.4-285416 


f 


3392 


< 5 


2952 




> 9 


174603.7-290247 


f 


1561a 


< 7 


3393 




> 3 


174606.2-290941 


gc 


3665 


< 11 


3392 




> 3 


174607.5-285951 


f 


2952 


■'*-16 


1561a 


209+2^^ 


6 


174610.8-290019 


gc 


3665 


< 4 


1561a 


""-13 


> 10 


174612.3-285706 


gc 


3665 


< 5 


0242 


78+" 

' °-15 


> 13 


174613.9-285924 


gc 


3393 


< 4 


2951 


47+25 


> 7 


174614.0-290220 


f 


3665 


< 5 


0242 


74t\l 


> 11 


174615.9-290257 


f 


2951 


< 11 


2943 




> 2 


174624.4-285712 


f 


3393 


25l? 


2954 


88iS 


3 



Table 8 

Hard X-ray Sources with Unusual Properties 





Fx" 




Continuum 




Iron Emission 


Periodicity 






Object 






r 


l-T 


Ff.. 


FM-'f,, 


P 


A 


Rcf. 




cr^ cm " ,^ ' ) 


(1(1-- cm 




(kc\ ) 


(kc\M 


lc\-) 




('<) 






AX .12315-0592 


50 


0.07 




17 


6.8 


900 


53()0 


90'= 




[1] 


RX J1802. 1+1804 


0.5 


13= 




> 7 


6.7 


4000 


6840 


lOO*" 




2 




AX J 1842. 8-0423 


4 


5 


2.9 


5.1 


6.7 


4000 








3 




XMM J174457-2850.3 


6.5 


6 


0.98 




6.7 


180 








'4' 




ICR J16358-4726 


70 


33 


0.5 




6.4 


130 


5580 


37 


[5,6] 


ICR J16318-4848 


8 


196 


1.6 




6.4 


2000 






[7,81 


ICR J16320-4751 


400 


21 


2.5 












[£ 


,10] 


XMM J174544-2913.0 


6.5 


12 






6.7 


2000 








[4] 


AX J1820. 5-1434 


23 


9.8 


0.9 




6.4 


90 


152 


57 




11 




AX J170006-4157 


5 


6 


0.2 






<1200 


715 


50 




12 




AX J1740. 1-2847 


4 


2.5 


0.7 






< 500 


729 


100 




13 




ISAX J1452. 8-5949 


0.6 


1.9 


1.4 




6.4 


> 1300 


437 


74 




14 




AX J183220-0840 


11 


1.3 


0.8 




6.7 


450 


1549 


63 




15 





"Observed flux, 2-10 keV. 

'^Modulation is only present below 2 keV. 

"^Partial-covering, intrinsic absorption. 
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